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Batch Studies for Methylene Blue Removal and 


Recovery by Untreated Coffee Residues 


Odysseas Kopsidas 


Laboratory of Simulation of Industrial Processes, Department of Industrial Management and Technology, University of Piraeus, 


Piraeus 18534, Greece 


Abstract: This paper reports the practicability of using coffee residues to remove widely used basic dyes like MB (Methylene Blue) 


from wastewaters. The effect of different batch system parameters, namely stirring speed, pH, initial dye concentration and contact 


time were studied. Moreover, in continuous fixed-bed column systems, the effects of parameters such as bed-depth, flow rate and 


initial dye concentration were examined. The experimental batch systems data were simulated using (a) Freundlich, Langmuir and 


Sips isotherm models and (b) Ist order, 2nd order, and intra-particle kinetic models. The results revealed that the MB is fairly 


adsorbed on coffee residues. This process could be a low cost technique for the removal of basic dyes from aqueous systems. 


Key words: Adsorption, desorption, column, MB (Methylene Blue), coffee residues, wastewaters. 


1. Introduction 


Many industries, such as paper, plastics, food, 
printing, leather, cosmetics and textile, use dyes in 
order to color their products [1]. In textile industries 
about 10-15% of the dye gets lost in the effluent 
during the dyeing processes which are harmful 
products and may cause cancer epidemics [2]. Dyes 
usually have a synthetic origin and complex aromatic 
molecular structures which make them more stable 
and more difficult to biodegrade [3]. The industrial 
runoffs are usually discarded into rivers and lakes, 
altering the biological stability of surrounding 
ecosystems [4]. Therefore, removal of dyestuffs from 
wastewater has received considerable attention over 
the past decades [5]. 

In wastewater treatment, various methods applied to 
remove dyes include photocatalytic degradation [6], 
membrane separation, chemical oxidations [7] and 
electrochemical process. Among the above mentioned 
techniques of dye removal, the process of adsorption 


gives the best results as it can be used to remove 


Corresponding author: Odysseas Kopsidas, Ph.D., main 
research field: environmental economics. 


different types of coloring materials [8]. 

Adsorption onto activated carbon is the most 
widespread technology for the removal of pollutants 
from water and wastewaters. The disadvantage of 
activated carbon is its high cost [9]. Hence, it is of 
pivotal importance thence of low-cost substitute 
absorbents to replace activated carbons. Various types 
of untreated biomass have been reported to have a use 
in dye removal: peanut husk, palm kernel fibre [10], 
Turbinaria turbinate alga, graphene, defatted jojoba 
and sugar beet pulp [11]. The sorption efficiency of 
several adsorbents is presented in Table 1. 

Further, 


materials are used to remove dyes in water and 


numerous pretreated lignocellulosic 
wastewater. Pyrolyzed date pits, date stones and 
Turbinaria turbinate alga have proved to be effective 


adsorbents. 


2. Materials 


Coffee residues and others similar byproducts like 
husks and 


investigated to remove various dyes and heavy metals 


coffee coffee grounds have been 


from aqueous solutions. Coffee husks have been used 
for the removal of Cu (ID, Cr (VD, Cd (IN and Zn (II), 
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coffee grounds for Nylosan Red N-2RBL and 
degreased coffee bean for Malachite Green. It is 
reported that the removal of Remazol Brilliant Blue 
RN and Basic Blue 3G is using coffee residues, 
focusing on batch processes and not on continuous 
systems [12]. 

In the current study, the adsorption capacities of 
coffee residues to remove MB (Methylene Blue) from 
aqueous solutions were investigated. Both batch and 
continuous | fixed-bed-column 


systems were 


investigated. As regards the fixed-bed-column 
systems filed with coffee residues, various bed-depths, 
flow rates and initial dye concentrations were 


examined for the first time. 


3. Methods 


Coffee residues, a low-cost material acquired from 
different cafeterias in the city of Piraeus (Greece) 
were used as adsorbent. This adsorbent was washed 
with distilled water and dried at 110 °C for 24 hrs to 
remove the humidity. 

The dye used herein in batch and column 
experiments was MB (Cj6H;sCIN3S°3H2O, molecular 
weight = 373.90 x 10° kg-mol') supplied by 
Sigma-Aldrich. A stock solution was prepared by 
dissolving a specific amount of MB (humidity 22%) 
m distilled water. Working solutions were 3-140 
mg'L!. MB concentrations were determined by 
measuring the absorbent values in each experiment 
with HACH DR4000U UV-VIS spectrophotometer at 
A = 664 nm. 

Isotherms were obtained from batch experiments. 
Accurately weighted quantities of approximately 0.5 g 
coffee residues were transferred to 0.8 L bottles, 
where 0.5 L adsorbate solution were added. The 
temperature was 23 °C, the MB initial concentration 
ranged from Cy = 1.5 mg'L” to 150 mg' L". The 
bottles are sealed and mechanically agitated for a 
period of 7 days. The 7-day period was determined 
after optimization analysis, with agitation periods 


varying from 4 hrs to 14 days, to ensure that nearly 


equilibrium conditions were achieved. The final 
concentrations were determined. Concentrations 
(before and after equilibrium) from each bottle 
represented one point on the adsorption isotherm plots. 
The pH was near 7.5. 

Batch experiments were carried out at stirring 
speeds from 65 to 664 rpm, pH varied from 1.54 to 
12.9, the adsorption temperature varied from 23 to 
70 °C, initial dye concentration from 3 to 140 mg'L” 
for a maximum contact time up to 95 mim. Samples 
were taken at 5 minutes intervals and the MB 
concentration was measured. The reactor, containing 
V= 1 L aqueous solution of dye was placed in a water 
bath to maintain constant temperature at the desired 
level. 

The study of the coffee residues by SEM (Scanning 
Electron Microscopy) was conducted at the Institute 
of Materials Science, National Centre for Scientific 
Research ‘Demokritos’, using an FEI Inspect SEM. 

The concentration of output solution was 
measurement at 4 = 664 nm and using HACH 
DR4000U UV-vis spectrophotometer. Finally, pH 
measurements were made using a digital pH meter, 


Multi Lab model 540. 


4. Results 


The SEM micrographs for coffee residues before 
and after MB adsorption are presented in Fig. 1. The 
coffee residues particle in Fig. 1(a) with magnification 
7,500X shows some pores. The surface texture of this 
particle 1s relatively rough as it can be observed in Fig 
l(b) with magnification 30,000X. The adsorption 
conditions were Co= 140 mg'L”, t = 95 min, mV = 1 
gL”, agitation speed = 664 rpm. The texture of the 
coffee residues particles after MB adsorption in Fig. 
1(d) is rougher comparing to the same material before 
MB adsorption, presented in Fig. 1(b). This fact 
indicates the swelling effect on the lignocellulosic 
particles after MB adsorption. 

Three isotherm models including Freundlich (1906), 
Langmuir (1916) and Sips (1948) [13] equations were 
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Fig. 1 SEM analysis for coffee residues. 


used to fit the experimental data. Freundlich isotherm 
model assumes that the surface of absorbent is where 
heterogeneous and polymolecular layer adsorption 


takes place. This model can be described byEq. (1): 
1 

q=K,-(C,)" (1) 
where g is the amount adsorbed per unit mass of the 
adsorbent (mg-g’'), C, is the equilibrium concentration 
of MB adsorbed (mg-L"'), Kr [(mg-g')(L-mg”)""] is 
the Freundlich isotherm constant related to adsorption 
capacity and n is the Freundlich isotherm constants 
related to adsorption intensity. In cases where the 
isotherm experimental data approximates the 
Freundlich equation, the parameters Kr and n can be 


estimated either by plotting logg versus logC, either 





j 
+ 


9/19/2012 HV | mag WD HFW | tilt 
12:36:28 PM | 10.00 kV| 30 000 x [10.0 mm|9.01 um} 0 ° 


9/19/2012 | HV | mag | WD | HFW EE 
1:00:30 PM |10.00 kV | 30 000 x |10.3 mm|9.01 uml 0º 


(d) 


by NLRA (Non-Linear Regression Analysis). The 
linear form of the Freundlich isotherm model can be 
defined by Eq. (2): 


l 
logg=logKp +—logC, (2) 
n 
The Langmuir (1916) isotherm model is given as: 
K 
q = LImCe (3) 
|+ KçC, 
Moreover, this equation in linearized form is: 
l l l l 
-=| —|+ — >| 8 
q \Im) (Ki Im) Ce 


where Kz is the Langmuir constant related to the energy 


of adsorption (L‘mg”) and qm is the amount of MB 
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adsorbed (mgg) when saturation is attained. In cases 
where the isotherm experimental data approximates the 
Langmuir equation, the parameters Kz and qm can be 
estimated either by plotting 1/g versus 1/C, either by 
NLRA. 

Numerous studies incorporate another important 
parameter, Rz, namely the separation factor. The value 
of Rz indicates the type of the isotherm to be either 
unfavorable (Rz > 1), linear (Rz = 1), favorable (0 < Rz 
< 1) or irreversible (Rz = 0) and is expressed by Eq. 
(5): 

l 


= —— — 5 
“ 14K,-C, e 


where Co is the initial dye concentration (mg: L’) and 
K, is the Langmuir constant (L:mg’). 

Sips (1948) [13] isotherm model is a combination of 
the Langmuir and Freundlich isotherm type models and 
is expected to describe heterogeneous surface better. 


The Sips equation is presented by Eq. (6): 


= Im (Ky C" 


6 
ak ae) o 


where Kz (L-mg 9, Am (mg-g "is the Langmuir 
constants and n the Freundlich constant. These 
parameters can be estimated by NLRA. 

MB adsorption isotherms experimental data for 
untreated coffee residues are presented in Fig. 2. The 
theoretical curves are estimated for the above 
mentioned three models. The parameter values of 
these models are shown in Table 1. These parameters 
were obtained by NLRA. The SEE (Standard Error of 


Estimate) was calculated in each case as Eq. (7): 


SEE = 





>, g Vaa) /(n'—p') (7) 
i=] 


where y; 1s the experimental value of the depended 
variable, Vitheor IS the theoretical or estimated value of 
the depended variable, n” is the number of the 
experimental measurements and p’ is the number of 
parameters (the difference n — p’ being the number of 
the degrees of freedom. 

The fitting of the Sips adsorption model to the 
present data was the most satisfactory for MB 
adsorption, better than the other two isotherm models, 
as shown by the corresponding SEE-values given in 
Table 1. The fitting of the Langmuir’s adsorption 
model to the present data was also satisfactory but to a 
lesser degree than the Sips model. The Rz values were 
found to be 0.236-0.935, i.e. 0 < Rz; < 1 for all MB 
concentrations Cy in the range of 3-140 mg-L". 

In order to make adsorption mechanism more clear, 
the adsorption data were fitted using the pseudo-first 
order, the pseudo-second order and the intra-particle 
diffusion kinetic model. 

The Lagergren (1898) [12] non-linear pseudo-first 


order equation is given as Eq. (8): 


q-g=q:e"" (8) 


where q and q, (mg: g') are the amounts of MB dye 
adsorbed per unit mass of the adsorbent at equilibrium 
time (f— œ) and adsorption time f¢ (min), 
respectively, while k (min is the pseudo-first order 


rate constant for the adsorption process. Moreover: 
q=(C-C)V/mand g,=(QG-CWV/m (9) 


where C, Co and C, (mg-L”) are the concentrations of 
MB in the bulk solution at time ¢, O and ©, 
respectively, while m (g) is the weight of the adsorbent 
used, and V (mL) is the solution volume. Further 


modification of Eq. (9) in logarithmic form gives: 


Table 1 Isotherms’ parameters of Methylene Blue dsorption on coffee residues. 


K-[(mg-g')(L:mg")'"" n 
Freundlich 3.555 1.612 
Langmuir 


Sips 1.195 


K; (L:mg"') dm (mgg"') SEE 
2.397 

0.0230 78.90 2.197 

0.0111 109.27 2.168 
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Fig. 2 Isotherm’s parameters of adsorption of coffee 
residues. 


In(q—q,)=Ing—k-t (10) 


The pseudo-second order equilibrium adsorption 


model equation is given as Eq. (11): 


E l ] 
q=q-|g +kot| or 4 =4-7 (11) 
— + kzt 
q 


where kz (min!) is the rate constant of second order 
adsorption. 

The intra-particle diffusion model based on the 
theory proposed by Weber and Morris (1963) [9] is 
expressed as Eq. (12): 


q=ct+ky, vt (12) 


where k, [mg (gminy! is the intra-particle diffusion 
rate constant and c (mg-g’ )is a constant related to the 
thickness of boundary. A value of c close to zero 
indicates that diffusion is the only controlling step of 
the adsorption process. 

The stirring speed was studied between 65 and 664 
rpm at constant dye concentration 14 mg-L", pH 
equals 7 and stirring time of 95 min. The correlation 
between the first-order rate constant k for the MB 
absorption by coffee residues is shown in Fig. 3. The 
rate constant is generally increasing with stirring 
speed up to 200 rpm were a level-off is reached. 

According to the literature, the pH of the aqueous 
solution is one of the most important parameters of the 


adsorption process. The effect of the initial pH on the 


MB adsorption onto coffee residues is illustrated in 
Fig. 4. The range of pH studied was from 2 to 12. As 
can be seen the rate constant increases linearly by pH 
increasing. 

The effect of initial MB concentration on adsorption 
on coffee residues is shown in Fig. 5. The kinetic 
parameters of the above mentioned three kinetic 
models were estimated using NLRA and are presented 
in Table 2. As can be seen, the SEE values are lower for 
the intra-particle diffusion kinetic model indicating 
better fitting to the experimental data. In addition the 
pseudo-second order kinetic model had better fitting 
that the pseudo-first one. 

The temperature variation from 20 °C to 70 °C on 
MB uptake by coffee residues revealed that increase in 
temperature of the process enhanced better adsorption 
of MB from bulk solution. The conditions under which 
they became the experiments with different temperature 


were pH = 8, the initial concentration was 14 mg:L", 


0,031 - 
0,029 - . 
0,027 - a 

— 0,025 | . 

= 0,023 - 

x 0,021 | 
0,019 - 
0,017 - 
0,015 | | | | 
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Fig.3 Effect of stirring speed. 
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Fig. 4 Effect of pH. 


492 Batch Studies for Methylene Blue Removal and Recovery by Untreated Coffee Residues 





+ 140 mg/L 
o 75.3 mg/L 
a 26.3 mg/L 
o 11.3 mg/L 
= 7.8 mg/L 
e 3.0 mg/L 




















t (min) 


Fig. 5 Effect of initial MB concentration. 


the stirring speed was 600 rpm and the dose of coffee 
residues was 1 g. Increase in temperature increased 
entropy of the system which yielded more chances of 
adsorption. The change in adsorption rate constant k as 
a function of adsorption temperatures is presented in 
Fig. 6. 
The rate constant k of the Lagergren model follows 
the Arrhenius law: 
k= p-exp(—E/ RT) (13) 
The activation energy E (kJ-mol’) for the adsorption 
of MB on untreated coffee residues can be estimated by 
linear regression of /nk on 1/T. This activation energy 
was found to be 27.4 kJ-mol” or 6.57 kcal-mol'. The 
frequency factor of Eq. (13) was p = 1,236 min”. 
Consequently, the intra-particle diffusion, which is a 
physical process, is probably the controlling step of the 


adsorption process. The parameters of fixed—bed 


column systems for coffee residues are shown in Table 3. 
This study is a research on the adsorption capacity of 
the coffee residue and its application to industrial scale. 
There is object for further research to determine the 
coffee residues quantity by region distribution and 
corresponding exploitation by industrial scale unit. 
The MB 


lignocellulosic materials found in the literature was 


adsorption capacity for various 
compared to the coffee residues adsorption capacity 
estimated in the present work (Table 4). The 
adsorption capacity of the coffee residue was better 
than the adsorption capacity of other waste biomass 
such as peanut husk, turbinaria turbinate alga, 
and wheat straw. The adsorption capacities of 
various coffee residues/byproducts found in the 
literature are shown in Table 5. The coffee residues 
adsorption capacity found in the present work was 
close to the average of the Freundlich and Langmuir 
capacities, 1.e., Kr and qm respectively, reported in the 
literature. 

The present study is a part of continues involvement 
of our Research Group in the study of wastewater 
treatment using original and modified (pretreated) 
lignocellulosic biomass as potential adsorbents. 
Previous study of our Group gave Langmuir capacity 
Gm = 38.7 mg: g! for pine sawdust, significantly lower 
comparing to qm = 78.9 mg: g! for the coffee residues 


studied in the present work. 


Table 2 Kinetic parameters of Methylene Blue adsorption on coffee residues. 


Cy (mg:L"') 140 75.3 26.3 11.3 7.8 3.0 
First order 

k 0.0183 0.0227 0.0211 0.0113 0.0107 0.0176 
qe 28.78 22.26 11.11 5.79 4.98 1.82 
SEE 1.2827 0.9917 0.7150 0.1271 0.0515 0.0243 
Second order 

K, 0.000376 0.000604 0.001139 0.000796 0.000823 0.004587 
qe 37.81 30.99 15.42 9.27 8.18 2.74 
SEE 1.1818 0.8594 0.6563 0.1198 0.0488 0.0209 
Intra-particle diffusion 

c -1.05 -0.70 -0.29 -0.49 -0.49 -0.16 
k, 2.3371 2.1586 1.0365 0.4275 0.3647 0.1706 
SEE 0.9196 0.6165 0.5137 0.1945 0.1613 0.0535 
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Fig. 6 Effect of temperature. 
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Table 3 Column models parameters of Methylene Blue adsorption on coffee residues. 


Q (mL/min) 


x= 15cm 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 

x=25 cm 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 
Bohart-Adams 
Clark 


C; (mg/L) 


165 
165 
130 
130 
150 
150 
145 
145 
800 
800 
1,600 
1,600 
1,600 
1,600 


150 
150 
800 
800 
550 
550 


20 
20 
40 
40 
70 
70 
80 


n 


2.000 
1.612 
2.000 
1.612 
2.000 
1.612 
2.000 
1.612 
2.000 
1.612 
2.000 
1.612 
2.000 
1.612 


2.000 
1.612 
2.000 
1.612 
2.000 
1.612 


N (mg/L) 


42,798 
40,522 
29,317 
26,317 
31,482 
28,382 
24,745 
20,673 
46,166 
33,869 
30,224 
22,024 
36,641 
25,017 


38,290 
35,501 
21,788 
16,007 
32,847 
30,028 


K 


0.000080 
0.000060 
0.001400 
0.000100 
0.000250 
0.000220 
0.000210 
0.000180 
0.000020 
0.000010 
0.000010 
0.000010 
0.000020 
0.000020 


0.000160 
0.000130 
2.25E-05 
1.98E-05 
8.61E-05 
7.51E-05 


Table 4 Methylene Blue adsorption capacity for various lignocellulosic materials. 


Adsorbents 


Activated carbon 


fibers 
Cocoa 


(The obroma cacoa) 


shell 
Coffee residues 
Commercial 


activated carbon 


Cotton stalk 
Cotton stalk 


Pretreatment 


Chenmical vapor 
deposition modified 
carbonized, burn 


activated 


sulfuric acid treated 


Krl(mg'g )(Lemg5!” 


120.000 
40.850 


3.555 


50.440 
202.210 


n 


1.460 
8.090 


qm(mg' g’) 
478.00 
212.77 


78.90 
370.40 


147.06 
555.56 


Qo (mg/g) 


96.90 
91.76 
66.39 
59059 
71.29 
64.27 
56.03 
46.80 
104.54 
76.69 
68.44 
49.87 
82.97 
56.64 


86.65 
80.34 
49.30 
36.20 
74.30 
68.00 


SEE 


2.3110 
2.3053 
0.4876 
0.6402 
5.8235 
5.1515 
T9229 
6.3513 
38.2861 
32.5900 
102.8607 
93.8242 
103.0998 
903.4492 


4.0545 
3.7795 
62.5900 
56.0000 
17.6150 
14.6410 


K,(L-mg"') 
0.3750 
0.2690 


0.0230 
1.0500 


0.0249 
0.6207 
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Table 4 continued 


Adsorbents Pretreatment Kil(mg:g (Lmg 5!” 
Cotton stalk prespnonçaçia 157.910 
treated 

Date pits activated Pyrolyzed, FeCl, 

74.871 
carbon activated 
Date stones activated Pyrolyzed, ZnCl 

71.190 
carbon activated 
Enteromorph a 
prolifera activated ZnCl, activated 236.230 
carbon 
Graphene - 90.920 
Jojoba (defatted) - 1.132 
Palm kernel fiber - 8.670 
Peanut husk - 16.650 
Pine sawdust Autohydrolyzed 15.680 
Pistachio hull - 112.300 
Posidonia oceanic 
(L.) dead leaves Carbonized, Z0Ch 112.120 

, activated 

activated carbon 
Sugar beet pulp - 12.520 
Sugar beet pulp - 2.411 
Tea (rejected) NaOH modified 40.870 
ee - 14.000 
turbinaria alga 
Turbinaria 
turbinaria alga pyres ay 
Tı urbinaria pyrolyzed, physically 105.000 
turbinaria alga activated 
Turbinaria pyrolyzed, 
turbinaria alga chemically activated FO 
Wheat straw Acid hydrolyzed 10.950 


n 


16.150 


4.670 


3.440 


3.720 


5.710 
1.054 
1.770 
3.270 
2.582 
4.200 


2.740 


1.990 
10.128 
2.110 


2.310 


3.030 


3.630 


5.670 
1.010 


Table 5 Adsorption capacity of various coffee residues/byproducts. 


Adsorbents 
Coffee grounds 


chemically activated 


carbon 
Coffee grounds 


chemically activated 


carbon 


Coffee husk-based 


activated carbon 
Coffee husks 


Coffee grounds 
Coffee bean 
(degreased) 


Coffee residues 


Coffee residues 
Coffee husks 
Coffee husks 
Coffee husks 
Coffee husks 
Coffee residues 


Adsorbates 


Nylosan Red 
N-2RBL 


Methylene blue - 


Remazol brilliant 
Orange 3R dye 
Methylene blue 2.567 


Methylene blue 51.926 
Methylene green 2.03 


Remazol brilliant 


Blue RN ne 
Basic Blue 3G 36.61 
Cu(II) 3.702 
Cd) 3.188 
Zn(II) 2:12] 
Cr(II) 2.353 


Methylene blue 3.555 


Krl(mg'g (L:mg’)" n 


qm(mg'g’) 


222.22 


249.46 


369.38 


270.27 


153.85 
167.00 
95.40 
72.13 
88.02 
389.00 


217.390 


244.600 
714.290 
242.110 


63.000 


163.000 


411.000 


345.000 
20.410 


dm(mg'g'') 


367 


181.8 


66.76 


90.09 
18.73 


55.3 


175 


240 

7.496 
6.854 
5.565 
6.961 
78.90 


Ku(L-mg”) 


14.516 


0.0830 


0.0429 


1.2300 


1.4400 
0.0009 
0.0317 
0.0850 
0.1072 
0.0960 


1.7000 


0.0149 
0.0039 
0.1300 


0.1830 


0.1040 


0.0830 


0.136 
0.663 


Kı(L:mg’) 


0.077 


0.109 


0.401 


0.0226 
0.2687 


0.0935 


0.009 


00.041 
0.4232 
0.3124 
0.2238 
0.1805 
0.0230 
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5. Conclusions 


The use of coffee residues is an innovative technique 
using waste biomass from urban and rural areas in and 
Industrial Ecology framework. Also, that is an index of 
the sensitivity of citizens for the rational management 
of waste and their contribution to sustainable 
development. The present work proves the viability of 
using coffee residues to remove basic dyes like MB 
from wastewaters. The effect of different system 
parameters on (a) batch and (b) continuous fixed-bed 
column systems were studied. The experimental 
systems data were simulated using the most commonly 
used isotherm and kinetic models. The continuous 
fixed-bed column results proved that the MB is 
practically adsorbed on coffee residues giving 
maximum Bohart and Adams capacity N = 46,166 
mg'L” or go= 104.5 mg: g! for bed-depth 15 cm, initial 
dye concentration 800 mg'L” and flow rate 20 
These 


suggesting this low cost residue for the removal of 


| . y 
mL'min . results provide evidence for 


basic dyes in industrial scale applications. 
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Abstract: This study sought to establish the compatibility effects of biochar/vermin-compost application on growth and yield of 


maize and cabbages. Biochar application rate was at 3 to 4% of soil weight. Biochar was mixed in various proportions to come up 


with treatments which were 100% vermicompost, 50/50% biochar: vermicompost, 75/25% biochar: vermicompost and 100% biochar. 


The experiment was repeated twice and arranged as a completely randomized design and replicated five times. The study established 


that maize yield increased with application of biochar from 25% to 50%. In the cabbage trial, comparable results in plant height, leaf 


numbers and final yield was obtained with 100% vermicompost application, 50% and 75% biochar inclusion. When biochar was 


applied at 100%, noticeable reduction in plant performance was noted. It is recommended to use biochar as a bio-fertiliser at 50/50% 


biochar: vermicompost in maize and cabbage production. 


Key words: Biochar, vermicompost, crop productivity, compatibility. 


1. Introduction 


In most arid and semi-arid regions, large areas of 
cultivated soils are characterised by high sand and low 
organic matter contents which in turn negatively 
affects soil physico-chemical properties [1]. These 
problems are particularly intense in developing 
countries like Zimbabwe where more than 80% of the 
population of smallholder farmers still rely on 
traditional technologies and tools; mainly hand-held 
hoes, minimal use of animal traction and no tractors 
[2]. However, soil health is the foundation of vigorous 
and sustainable crop production. This is characteristic 
of the semi-arid tropics like Zimbabwe. As a potential 
solution, the use of biochar as a soil amendment has 
been introduced. 

Biochar is a carbon rich material produced by 
different conversion 


thermo-chemical biomass 
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process like pyrolysis. Pyrolysis 1s a dry carbonization 
technique typically operated at 400-800 °C [3]. 
The application of biochar to the soil has been shown 
to enhance the capacity of the soil to retain 
nutrients and lessen leaching [4]. Studies done by 
Githinji, L. [5] have revealed that addition of biochar 
increases soil porosity and aggregation. Biederman, L. 
and Harpole, W. S. [6] reported that biochar addition 
to the soil leads to increased above ground 
productivity, crop yield, soil microbial biomass, 
rhizobia nodulation, plant potassium concentration 
and total soil nitrogen. Jeffery, S., et al. [7] found a 
10% increase in the mean yield of various crops. The 
percentage increases found 1n literature tends to differ 
and depends on various factors such as initial soil 
High 
applications of biochar have been reported to reduce 
yields [9]. 

A study by Carter, S., et al. [10] reported increases 


properties and biochar characteristics [8]. 


in above ground biomass after biochar additions to 


non fertilized soils. However, other studies have found 
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that biochar increases the persistence of simazine 
herbicide. Despite the promotion of biochar as a soil 
amendment, its compatibility with other soil 
amendments such as vermicompost 1s currently 
unknown. Vermicompost is one of the soil organic 
amendments that has found its way into Zimbabwean 
fields and gardens and now that biochar is coming in, 
the compatibility of these amendments need 
determination. The study explored the effects of 
biochar mixes on the 


various vermicompost 


productivity of cabbage and maize in a pot trial. 
2. Materials and Methods 
2.1 Experimental Site 


Two experiments were carried out at BUSE 
(Bindura University of Science Education) nursery 
and were repeated twice. BUSE is located at the 
coordinates: 17°18’58’’ S and 31°19°23” E. The two 
sites are located in agro-ecological region lla of 
Zimbabwe [11] and receive about 800 mm rainfall per 
annum. The average temperature for both areas was 
25 °C. The soil used at both sites was medium grained 
sands with a pH of 4.1 (CaCl,). 


2.2 Experimental Details 


Biochar used in the experiments were obtained 


from the pyrolysis of maize stova. 


Biochar/vermicompost mixes used in these 
experiments were in the proportions 50/50, 75/25 
vermicompost/biochar and 100% biochar. The mixes 
were used to fill the pots. Maize fert (7% N, 14 P, 7 K) 
was applied at a rate of 300 Kg/ha. At topdressing, 5 
grammes of ammonium nitrate (35.5%) was applied at 
a rate equivalent to 250 kg/ha. The experiment was 
arranged as a completely randomized design with five 
replications. The experiment was repeated twice over 
time. Watering was done at 50% moisture depletion 
calculated using the evaporation pan measurements 
relevant at each crop stage. For tomatoes the weekly 
fungicide sprays were done through alternating copper 


oxychloride and Dithane M45. 


2.3 Data Collection and Analysis 


The data from the two repeats of the same 
experiment were pooled together prior to analysis. The 


data was analysed according to the model: 


Y. =u+T +e, 


y 


Where u is the average mean, T; means the 
treatment effects and e; is the error term. 

Data on plant growth parameters were subjected to 
analysis of variance using Genstat Version 13. The 
(LSD (Least 


Difference) 0.05) was used to separate treatment 


least significant test Significance 
means of the measured parameters with significant 


differences. 
3. Results and Discussion 


3.1 Maize Leaf Number 


The results indicate that biochar inclusion at 50% 
performed similarly to vermicompost but was 
significantly (p < 0.05) higher compared to inclusion 


at 25% (Fig. 1). 
3.2 Maize Yield 


The same trend was observed on maize yield with 
biochar inclusion at 50% having the highest yield 
compared to the rest of the treatments (Fig. 2). 


3.3 Cabbage Height 


There was no significant variation in cabbage plant 
height when biochar is applied at 50% and 75% when 
compared to vermicompost alone (p > 0.05). The 
application of biochar in cabbage production at 100% 
resulted in significantly reduced plant height (p < 0.05) 
when compared with application of vermicompost 
alone or biochar at 50% and 75% (Fig. 3). 


3.4 Cabbage Leaf Canopy 


The cabbage leaf canopy diameter was similar at 50% 
and 75% biochar inclusion when compared with 100% 
vermicompost application (p > 0.05). However, at 


100% biochar application rate, there was a significant 
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Fig. 1 Effect of various biochar/vermicompost mixes on maize leaf number. 
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Fig. 2 Effect of biochar vermicompost mixes on maize yield. 
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Fig. 3 Effect of biochar vermicompost mixes on Cabbage plant height. 
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Fig.4 Effect of biochar/vermicompost mixes on cabbage leaf canopy. 


Cabbageb Leaf 
number 





Fig.5 Effects of biochar/vermicompost mixes on cabbage leaf number. 


reduction in the cabbage leaf canopy (p < 0.05) 
(Fig. 4) 
3.5 Cabbage Leaf Number 


There was no significant difference in cabbage leaf 
numbers grown under vermicompost alone compared 
to 50% and 75% biochar application rate (p > 0.05). A 
significant reduction in the number of cabbage leaves 
was noted when biochar was applied at 100% 
compared to vermicompost only and biochar 


application at 50% and 75% (p < 0.05) (Fig. 5). 


4. Discussions 


The study focused at establishing the effect of 


Biochar/Vermin-compost application on growth and 
yield of maize and cabbages. Maize yield increased 
with increasing application of biochar from 25% up to 
50%. This is in tandem with the results of other 
researches from elsewhere in the world [12, 13]. The 
noticeable higher performance in maize with increased 
biochar application could be attributed to the 
improved soil physico-chemical characteristics of the 
soil which include among other properties the capacity 
of the soil to retain nutrients, increased cation 
exchange capacity and lessened leaching. 
Vermicompost could be the supplier of the nutrients 
while biochar preserves them for crop usage. The 


results also serve to confirm the compatibility of the 
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two organic soil amendments on crop productions. 
The usage of the mixes may be limited to the 50% 
inclusion level as any further increase in biochar did 
not yield a proportional benefit. 

This agrees with Yamato, M., et al. [14] who 
indicated that the mobile portion consists of biochar 
short-term effects, the resident portion carries its 
long-term effects which vary with concentration levels. 
Within the mobile portion, the organic part acts much 
like other degradable carbon sources (e.g. compost or 
detritus). It contains dissolved organic carbon and 
organic matter that is available to soil microbes. The 
availability of the organic matter for plant growth will 
depend on the quantities applied to the soil. 

In the cabbages, there was a similar response 1n all 
parameters measured to vermicompost/biochar 
application, comparable results in plant height, leaf 
numbers and final yield were obtained with 100% 
vermicompost application, 50% and 75% biochar 
inclusion. When biochar was applied at 100%, 
noticeable reduction in plant performance was noted. 
This agrees with Verhejien, F., et al. [15] and 
Lehmann, J., et al. [16] who established that biochar 
helps in nutrient retention and cannot work on its own 
in nutrient supply. Biochar works in combination with 
other soil components such as soil microbes to 
improve the overall soil dynamics, the conditions in 
which plants grow over a long term. This allows for 
better plant nutrition, improved plant growth and yield 
and brings overarching benefits to agricultural soil 
productivity. 

This was further confirmed by Vincent, V., Thomas, 
R. G., Chan, K. Y., Xu, Z. and Julie, M. [11, 17, 18] 
who noted that, one reason for the ability of 
Amazonian Terra Preta soils to maintain high fertility, 
compared to adjacent infertile soils, is their ability to 
retain nutrients. 

Biochar has properties for high cation exchange 
capacity which allows it to retain positive ions in a 
form that is exchangeable. These positive ions, vital 


nutrients such as ammonium or potassium cations, are 


availed to plants when needed. This is so, because 
oftentimes when conventional fertilizers are applied to 
poor soils, the nutrients end up lost through leaching, 
weeds, or microbial activity before crops utilizes 


them. 


5. Conclusion 


It can be concluded from this study that the there 


are synergistic effects of mixing biochar and 


vermicompost with a limit of 50% inclusion rate with 
vermicompost in order to improve maize and cabbage 


productivity. 
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Abstract: The cave-site of Gruta da Furninha 1s a coastal karstic cavity. In the late 19th century, excavation of the sedimentary infill 
of the cave allowed the identification of two lithostratigraphic units: a Holocene one dominated by sands, containing human bones 
and a Neolithic industry; and a Late Pleistocene one, containing a diverse set of fossils set and Paleolithic artifacts. This study mainly 
focuses on taphonomic and palaeoenvironmental data concerning the Pleistocene bird bones that were collected from six layers (at 11 
m and 6.7 m asl). A main result of this study was the discovery of a left humerus fragment, with osteological features of the Alcidae 
family; from comparison with upper arms of distinct species of this family, it was concluded that this fossil belongs to Penguin 
impennis. The Pleistocene birds of the Furninha cave were also compared with the avifauna that currently occurs in the region of 
Peniche peninsula and a climate and environmental interpretation of the Pleistocene fossiliferous set is provided. 


Key words: Birds, Pleistocene, marine caves, Portugal. 


1. Introduction Zbyszewski, G. [2] identified at the base of the lower 


a unit (at 6 m asl) a marine stratigraphic level, which they 
The Furninha cave is a coastal karstic cavity in 

no ascribed to the Riss-Wurm interglacial. 
Lower Jurassic limestones, located ca. 850 m SE of 

In the upper unit “Entulho Superior” (Upper Rubble), 
Cape Carvoeiro, on the southern coast of the Peniche 
some of the bones collected are human. In the ca. 8 
peninsula (western coast of central mainland Portugal), 


about 75 km north of Lisbon (39°21°23” lat. N and 
9º26'14”” long. W). This cave has an entrance of 


marine genesis located 1n the middle part of the coastal 


m-thick lower unit “Areias Quaternárias” (Quaternary 
Sands), with 7 mam “níveis ossiferos” (bone-rich 
levels) I to VII, were recognized, several separated by a 
travertine and nine of them rich in well-preserved 
cliff, at ca. 15 m asl. co. 
faunal bones [1]. A description of the sedimentary 
The sedimentary infilling of the cave was excavated 
in the late 19th century by Delgado, J. F. N. [1], which 


defined two main stratigraphic units: a Holocene unit 


succession, including fossils and artifacts, can also be 
seen in Ref. [3]. 
Cunha, P. P., et al. [4] made a detailed 
and a Pleistocene unit. Later, Breuil, H. and o 
characterization of this coastal area. Several wave-cut 
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24-27 m; 20-21 m; 15-17 m; 10-12 m and 6-8 m. 

The first vertebrate inventory of fossils collected 
from the Pleistocene unit listed ca. three dozen 
mammalian taxa and showed the presence of birds, fish 
and chelonians [1]. Later, this Late Pleistocene fauna 
was more thoroughly studied by Harlé, E. [5]. In the 
case of birds, Harlé benefitted from the collaboration 
of Newton. Some other publications addressed the 
fossil remains of Furninha cave [6-8], but all of them 
just repeated the data already published by Harle. 
During the 20th century and early 21st century, the 
mammals were revisited, and new analyses were 
provided [9, 10]. In the case of birds, until 2007, the 
bibliography only reproduced the results of the study 
made by Newton in Harlé [5]. Since 2008, there have 
been some new studies (paleontology, taxonomy, 
taphonomy and paleoenvironment) focused on fossil 
birds from Furninha cave [9-15]. 

Until now, the bibliography related to the Furninha 
cave avifauna had two gaps: (1) the downgrading of 
taxonomic nomenclature; and (11) the lack of detailed 
study of the majority of the bird remains (only 20% of 


the total was previously studied). 
2. Material and Methods 


The present study is focused on the analysis of all the 
bird remains from Gruta da Furninha (436 bones), 
currently housed at the Geological Museum of LNEG 
(Laboratorio Nacional de Energia e Geologia) in 
Lisbon. The analyzed avifauna remains of Furninha 
cave were complemented with literature search and 
through comparison with the reference osteological 


collections (Arqueosciences Laboratory of General 


Direction of Cultural Heritage; Achaeozoology and 
Paleontology Laboratory of Portuguese Center of 
Geo-History and Prehistory). The association of 
taxonomic data with the environment of the species 
identified, according to the respective stratigraphic 


level of provenance, was also carried out. 
3. Previous Studies 
3.1 Geoarchaeology 


Delgado, J. F. N. 


outstanding detail and the archaeological materials 


[1] published results in 


have been successively reviewed or acknowledged [2, 
3, 5-7, 16-19]. The upper lithostratigraphic unit 
(“Entulho Superior”) encompasses an early to late 
Neolithic necropolis with abundant organic material, 
osteological human remains, personal adornments, 
bone tools and pottery along with polished and 
knapped stone tools [18] (Table 1). 

The lower lithostratigraphic unit Quaternary Sands 
(“Areias 


succession 


Quaternárias”) comprises a sedimentary 
containing artifacts encompassing a 
sequence between the Middle and Upper Paleolithic, 
with assemblages pointing to the Acheulian, 
Mousterian, Gravettian, Solutrean and Magdalenian 
[3], most 
industries—in association with abundant and diverse 
faunal remains [1, 5-15] (Tables 1 and 2 and Fig. 1). 


Unfortunately, the site has very few absolute dates. 


of them—especially the Mousterian 


Two radiocarbon AMS (Accelerator Mass Spectrometry) 
dates, one on human bone and another on a bone pin 
from the upper unit (Holocene) confirm that part, at 
least, of the necropolis identified from the upper unit 
dates from the Late Neolithic [18, 20, 21]. Another age, 


Table 1 Stone-tool assemblage from Furninha cave, according to Delgado, J. F. N. [1]; Bicho, N. and Cardoso, J. L. [3]; 


Cardoso, J. L. and Carvalho, A. F. [18]. 


Flint Quartzite Quartz Limestone Amphibole Other Total Reference 
Acheulean l 1l 
Mousterian 25 2 4 31 3 
Gravettian 24 l 25 3 
Solutrean 9 9 3 
Magdalenian 8 8 3 
Indeterminate 46 l 47 3 
Neolithic 70 l l 26 4 102 3 
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Fig. 1 Stratigraphy, plan and localization of the Furninha cave. 


Table 2 Stratigraphic levels of the Pleistocene lithostratigraphic unit and correspondence of the different fauna groups for 


each of these levels, according to Delgado, J. F. N. [1]. 


Main subunits Vil VI V IV 
Stratigraphic 

Ei 1(13.5) 1013) 212.8) 1 00) 19.2) 
Mammals * * * * P 
Birds * * * 

Reptiles * * 

Amphibians * 

Invertebrates 

Artifacts * * x $ 


resulting from Uranium-series dating of a bone 
collected from the 3th stratigraphic level (8.5-9 m asl), 
is 80,880 + 42, 420 — 31, 260 years [22] (Fig. 1), but 
cannot be used as a reliable chronological marker due 


to the large standard deviation. 


3.2 Paleontology 


13-15] allowed the 


identification of the vertebrate taxa of the Furninha 


Previous studies [1, 5-10, 


cave faunas and their distribution by the various 
stratigraphic levels (Figs. 1 and 3): carnivorous 
mammals (Ursus arctos, Hyaenaprisca, Panthera 
pardus, Lynx pardina, Canis lupus, Vulpesvulpes, 
Melesmeles and Martessp); herbivores (cervids, 
Dicerorhinus hemitoechus, Bos primigenius, Equus 


caballus and Sus scrofa); Paleoloxodont antiquus, 


1 (9) 


x 


x 


m II I 
2(8.7) 3(85) 172) 201) 1() 265.7 
x * x * x * 

x * x * 
x 

* x 
* x * 


Oryctolagus, flying mammals (order Chiroptera); 


amphibians; turtles; lizards and birds. 


4. Results 
4.1 Osteological Data 


A total of 436 bones of birds were identified (10% of 
the total number of vertebrate remains), collected from 
six of the seven Pleistocene stratigraphic levels of the 
Furninha Cave (Table 2). Most (95.4%) are bones of 
the appendicular skeleton. Of these, 36.2% are from 
hind limbs and 59.2% are from anterior limbs (Fig. 2). 
Of the remaining bones, 4.3% belong to the axial 
skeleton and only 0.3% to the cranial skeleton (a skull 
of P. graculus). The remaining variables have opposite 
signs, indicating positive contribution for ozone 
creation [10, 12]. 
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59,2% 


Fig. 2 Quantity and percentage of the skeleton parts represented by bird bones of the Furninha cave. 


4.2 Taxonomy 


Although the species T. tadorna predominates, with 
149 bone specimens identified, Furninha cave has a 
large taxonomic diversity, when compared with other 
Portuguese Pleistocene sites. The predominant groups 
are the passerines and the anseriforms. Anseriforms are 
dominated by marine species (Table 3, Fig. 3), which is 


explained by the location of the cave. 
4.3 Paleoenvironmental Interpretation 


Except for the bone of P. impennis, which have a 
fracture that could be of anthropogenic origin, the bird 
bones from Furninha cave do not provide evidence of 
human activities (cut marks, effects of fire) [12, 14]. 
Damage from predation was also not identified. 
Although the species identified are not typical of 
caves, they are marine species and the Furninha cave 
is located on a sea cliff. It is believed that the 
accumulation of bird remains would have resulted 


from animals that died inside the cave, but also other 


remains were brought by other animals and by the 
primitive humans using this cave. So, the biological 
remains also contributed to the sedimentary filling of 
the cave, which mainly consists of sand brought by 
wind and gravels provided by gravity processes, 
which would explain the large proportion of broken 
bones. 

The Furninha cave avian assemblage is dominated 
by seabirds (anseriforms and pelecaniforms). However, 
the presence of continental species (such as the grey 
partridge) indicates a connection with continental 
environments. Species associated with the following 
environments were identified: coastal habitats, forests, 
open fields, woods, wetlands and mountain areas. 

The presence of birds as the P. impennis, T. tadorna, 
C. olor and S. mollissima point to a colder than 
present climate. Other bird species, such as A. crecca 
and N. phaeopus, as well as mammal species collected 
in stratigraphic association, suggest humid and hot to 


temperate climate conditions [9, 12, 13]. 
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Fig. 3 Distribution of major groups of mammals and birds by the different daterat: Pleistocene subunits of Furninha 
cave, according to Delgado, J. F. N. [1], Brugal, J. P., et al. [10] and Figueiredo, S. [13] 
On the right side, the letters indicate the type of climate conditions suggested by the bird species (Table 3). The blue line on the right 


represents the distribution of Paleolithic industries, indicated by the artifacts found in the various levels. The presence of 


Palaeoloxodon (a fragment of tooth lamella) was identified [13], but it does not have the indication of the layer of origin, reason why 


its representation in this figure was not placed. 


Table 3 List of bird species from Furninha cave, according to Brugal, J. P., et al [10] and Figueiredo, S. [12] 


Distribution of species by interpreted environment: A species can occur in more than one environment, in this case it was divided 1 


by the number of environments where the species occur (Bos = Woods; CA = Open field Areas; Zh = Wetlands; LAD = Freshwater 


lakes; Mont = Mountain areas; ZR = Rocky areas; Flo = Forest areas; Cos = Coastal Areas). Climate: A (wet-temperate), B 
(wet-cold), C (dry-temperate), D (dry-cold). 


Anseriforms 


Galiforms 


Columbiforms 


Charadriiforms 


Tadorna tadorna 


Tadorna ferruginea 


Somateria mollissima 


Somateria sp 
Cygnus olor 

Anas crecca 

Anas sp 

Melanita nigra 
Alectoris rufa 
Coturnix coturnix 
Columba livia 
Numenius phaeopus 
Pinguinus impennis 
Gallinago sp 


Larus sp 


CA Zh 
0.5 
l 
l 
0.33 
0.5 


Environments 
LAD Mont 
0.5 
1 


0.5 
0.5 


0.33 


0.5 


ZR 


0.5 


Cos 


0.33 


Climates 
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Table 3 continued 


Bos CA 
Procellariiforms Puffinus puffinus 
Pelicaniforms Phoenicopterus ruber 
Phoenicopteriforms Phalacrocorax aristotelis 
Tito alba 0.5 0.5 
Strigiforms Bubo bubo 
Ásio flammeus 0.33 0.33 
Falco tinuculus 
FOTO las chrysiateus 0.33 
Aquila sp 0.33 
Gyps fulvus 0.5 
Pyrrhoc. pyrrhocorax 
Pyrrhocorax graculus 
Corvus corone 0.33 0.33 
Corvus frugilegus 0.5 0.5 
Corvus monedula 0.33 0.33 
Passeriforms Turdus merula 0.5 0.5 
Turdus philomelos 0.5 0.5 
Turdus pilaris 0.5 0.5 
Turdus iliacus 0.5 0.5 
Turdus sp 0.5 0.5 
Pica pica 0.33 0.33 
Cuculiforms Cuculus canorus 0.33 0.33 
Average 0,42 0,48 


5. Discussion 


A relevant highlight is the discovery of a left 
humerus fragment found in the middle of other bones 
not classified. This bone has osteological features of 
the Alcidae family: pronounced longitudinal flat of the 
diaphysis, lateral and ventral projection of proximal 
humeral articulation and slightly prominent pectoral 
crest, compared to dorsal margin. From these features, 
and from comparison of upper limbs of distinct 
species of the Alcidae family, it was concluded that 
this specimen belongs to P. impennis [11]. 

Most of the Pleistocene birds identified belong to 
taxa still living in Portugal, but only the following 
species have ceased to occur in the peninsula and its 
surroundings: A. chrysaetus, C. frugilegus, G. fulvus, P. 
ruber, P. impennis, S. mollissima, P. perdix and P. 
graculus. There are still a few reports of S. mollissima 


from the Minho estuary (upper North coast) to the 


Environments 
Zh LAD Mont ZR Flo Cos 


Climates 


l B 

l A 

l A 

D 

0.5 0.5 D 

0.33 D 

0.5 0.5 C 

0.33 0.33 C 
0.33 0.33 C 
0.5 C 

0.5 0.5 D 

l D 

0.33 C 
C 

0.33 C 

D 

D 

D 

D 

D 

0.33 C 

0.33 C 


0,35 0,55 0,66 0,5 0,38 0,73 


Q 
x 


Algarve, although none in Peniche [23]. These reports 
must be checked by the PCRB (Portuguese Committee 
for Rare Birds). The same holds true for C. olor, which 
also occurs around Peniche, but most are probable 
escapees from captivity. 4. Chrysaetos individuals 
prefer high inland mountains, but have occasionally 
been spotted in the Serra d' Aire e Candeeiros, not far 
away. G. fulvus also prefers inland habitats, but there 
are records of stray individuals spotted in the Tagus 
estuary and in an urban environment near Oporto. C. 
frugilegus has also become a rare species in Portugal, 
only recorded occasionally in the inner south and in the 
Algarve. It is nonetheless a bird that is common in 
colder climates. P. graculus is completely absent from 
Portugal, but is very common in several Portuguese 
Pleistocene assemblages. It is a bird that lives in high 
mountains above the 1,300 m, which makes the 
presence of its bones in Furninha remarkable, not 


because it is close to the coast, but because of the 
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latitude and the altitude as well as the modern-day mild 
climate of the place. P. ruber, an American bird from 
the Caribbean region, has been separated from the 
European species, and is about the same size as the P. 
roseus, which in his turn is quite common in the Obidos 
lagoon. Although Svensson, L., et al. [24] tell of its 
presence in Europe, he also claims they are of captive 
origin; no record from Portugal is known. As older 
guides [25-27] make no distinction between them, it is 
assumed that the bones found represent P. roseus 
instead, which is much more plausible. 

The subunit HI (8.5-9 m asl) has the highest number 
of faunal remains (mammals and birds) and has the 
presence of species associated with different climates 
(warm: G. fulvus, A. rufa or D. hemitoechus and E. 
caballus; cold: P. impennis or P. graculus) (Fig. 3), 
suggesting that this stratum was formed over a period 
of transition from an interstadial to a colder period. 

The long period with successive faunal and human 
occupations of the cave occurred after the cave was 
above sea level, following the genesis of the marine 
platform at 15 m asl. The sedimentary filling of the 
lower lithostratigraphic unit (Pleistocene), which 
contains Middle Paleolithic to Upper Paleolithic 
industries (these only at the top) interspersed in aeolian 
sands, must encompass the period of ca. from 90 to 12 
ka. This is supported, although without great precision, 
by the only absolute date available. 


6. Conclusions 


Building a bridge to the present-day avifauna 
hereabouts, it can be concluded that it has barely 
changed since the Pleistocene. From the analysis of a 
satellite map from Google Earth, including the 
Berlengas archipelago and the surroundings up to 
Caldas da Rainha and down to the São Domingos Dam, 
and having visited the region, the diversity of 
ecosystems can be attested, both terrestrial and marine: 
open sea, islands, coasts with both sand beaches and 
high cliffs, the harbor itself, and on the mainland 


agriculture fields both big and small gardens, with a 


multitude of trees and bushes of several sizes, a few 
streams, the Obidos lagoon, the São Domingos Dam 
and the Paul da Tornada, a biodiversity hotspot. It 
comes as no surprise then that from a search on the 
eBird database for only two of the hotspots in the 
Peninsula of Peniche (the Carvoeiro Cape and the 
harbor itself, thus on both sides of the Furninha cave), a 
total of 162 different species have been identified from 
the most diverse groups. From the Atlas of Breeding 
Birds, in its already somewhat outdated second edition 
of 2008, a total of 119 species are reported to nidify on 
the area. Both lists include pelagic birds, coastal birds, 
waders, raptors, and all sorts of smaller birds that are 
common in Europe. 

The identification of P. impennis is important, 
because this is a species associated with a cold climate, 
and because this species is poorly known in the fossil 
record of Portugal: Until this identification, it was 
only known in Gruta da Figueira-Brava and in Porto 
Santo [28, 29]. 

The distribution of fauna (birds and mammals) in 
relation to sedimentary levels allows a study of 


corresponding paleoclimatic evolution [9, 10, 12, 13]. 
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Abstract: Microscopic analysis of the skin and cuticle of wool fibers in Dubska pramenka show different qualitative parameters. The 
research study included different parts of the body and the samples of the wool fibers from different regions of the body. Histological 
description shows difference in the basic structurers of the skin in Dubska pramenka. A very pronounced thin epidermis, while 
dermis and hypodermis are more developed on the samples from the shoulder. Cornified flakes-cuticle of the root of the tail was of a 
fine structure, and cornified flakes from the rump showed the features of the rough wool fibers. “Transitional form of cuticle” was 
dominant on the shoulder. The qualitative parameters have a significant influence on the overall quality of wool in general, and the 
research contributes to a greater usability value of the wool and development of livestock farming in areas in a broader sense. 


Key words: skin, cuticle, Dubska pramenka. 


1. Introduction studies that focused on the structure of the skin and 


o the characteristics of the cuticle found differences m 
Skin is the largest organ and it covers the entire 
i , , i the structure, especially in the outer layer of the 
body. It consists of epidermis, dermis and , , , , 
, , , skin—epidermis, and the appearance of the cuticle, 
subcutaneous tissue or subcutis. The structure of skin , , , 
, , , , i , which varied depending on the region of the body. In 
in domestic animals is not equal, primarily, the l 
, , , , , general, total quality of the wool depends on the 
epidermis-dermis ratio varies, the presence and o 
qualitative parameters of the skin and fibers and so 
thickness of subcutaneous connective tissue. : 
nas does its classification and theuse in the processing and 
Difference in histological structure of some areas of Na 
development of livestock farming 1n areas. 
the body is also present in the same animals [1]. Wool, 


which is a product of the skin, is made of a set of wool 2. Material and Methods 


fibers. The surface of the wool fibers contains co 
The samples for microscopic examinations of the 
cornified cells, flakes, i.e. cuticle. Katica, A., et al. [2] 
skin and fibers of Dubska pramenka were taken from 


Corresponde aline Nadia Niles. PhD, man Vlasic mountain in Central Bosnia. The animals were 


research field: histology in veterinary medicine. reared under nomadic conditions, at the latitude of 
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850-1,600 m. The animals were healthy with a good 
body condition, weighing between 50 and 65 kg, 
above three years of age. Nutrition of the sheep was 
usual for nomadic conditions during summer, with the 
addition of grains before estrus, in line with the 
research studies [3, 4]. The samples of the skin for 
histological analysis were taken from the middle of 
the last left rib, the shoulder, the rump and the root of 
the tail according to the Markotic method. A moment 
before sampling, the site was previously shaved, 
washed and disinfected, and skin cuts with subcutis 
were prepared, 6 mm in diameter. The samples were 
placed in plastic bottles with caps that contained 10% 
formalin until histological processing 1.e. molding into 
paraffin blocks. Histological preparations of the skin 
were prepared according to the usual sample 


preparation protocol for optical microscopy 
observations. The analysis of the microscopic skin 
preparations was performed using binocular optical 
microscope MOTIC TYPE 120M with magnification 
of 100, 200 and 400 times. The fiber for microscopic 
analysis was taken from different parts of the body: 
the shoulder, the rump and the root of the tail, cutting 
the wisps next to the skin. Prior to microscopic 
examination of the cuticle, all mechanical dirt is 
cleaned; the samples were washed with neutral soap 
and rinsed with water, and then rinsed with distilled 
water. After that, the samples were immersed in 
hydrogen peroxide for 24 hours, then in xylol for 48 
hours. These treated samples were rinsed with 
distilled water and left to dry. The processed samples 
were placed on the glass slide, embedded into glycerin 
and examined under the microscope. The fibers were 
analyzed with binocular optical microscope, with 


magnification of 100, 200 and 400 times. 


3. Results and Discussion 


Microscopic analysis of the Dubska pramenka skin 


was performed using epidermis, dermis and 


hypodermis. Using descriptive presentation, the 


particulars of the skin of Dubska pramenka were 


singled out from different parts of the body. 
Microscopic examination of the skin of the Dubska 
pramenka taken from the shoulder identified a very 
thin epidermis, with dermis and hypodermis well 
developed. The attention was focused on dermis 
and in that particular segment, that the presence of 
hair follicles grouped in sets of two or three follicles 
(one central and two peripheral) and rarely asingle 
follicle was found [5]. 


collagen and elastic fibers that give dermis the 


Reticular dermis contains 
appearance of compactness (Fig. 1). Sebaceous 
glands (Glandulae alveolares) are immediately next to 
hair follicles. Sweat glands are surrounded by the 
elements of connective tissue where adipocytes can be 
seen [6]. 

The analysis of the sample fibers from the 
shoulders shows a so called transitional form of 
cuticle (Fig. 2). The cornified flakes are rather large 
with multi-corners, and the clearly pronounced edges 
are merging forming a whole. 

A thin epidermis is characteristic on the rump of 
Dubska pramenka. Dermis is well developed, with 
connective tissues and notable nuclei of the smooth 
muscle cells. The compacted collagen and elastic fibers 
are present in the papillary layer, while hair follicles 
are surrounded with poorly developed sebaceous 


glands, which are globular or elongated (Fig. 3). 





Fig. 1 The skin of the shoulder (Dubska pramenka); the 
white arrow indicates dermis; the blue arrows indicate 
sweat glands; the green arrow indicates sebaceous gland; 
HXE, x 200. 
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Fig. 2 Native preparation; x 400; shoulder. 





Fig. 3 The skin of the rump (Dubska pramenka); the 
white arrow indicates epidermis; the blue arrow indicates 


hair follicle; the green arrow indicates sebaceous gland; 
HXE; x 200. 





Fig. 4 Native preparation; x 400; the rumps. 


On the wool fiber samples taken from the rump, 
cornified flakes overlap continuously (Fig. 4). They 
have irregular shape reminding of roof tiles. This 
appearance of the cuticle most frequently 1s associated 


with rough wool fibers. 


Epidermis of the root of the tail is notable, but 
without clear strata. Its surface shows a layer of 
squamous cells, darker in color, flat, interconnected. 
This layer makes the stratum corneum, made of cells 
without nuclei, with cytoplasmic inclusions of 
keratin. The skin of the root of the tail shows clearly 
visible outer cells separated from other layers that do 
not have clear edges, and they make stratum 
disjunctum. Dermis is rich with oval hair follicles 
with developed sebaceous glands. Loose connective 
tissue is present in the outer papillary layer, while 
deeper dermis layers also contain collagen fibers 
(Fig. 5). 

Cross sections of sebaceous glands are noted 
(Glandulae glomiformes), which is highlighted in the 
studies by Janquera L. C. [7]. The studies to date 
indicate that the distribution of sebaceous glands in 
this part of the skin may affect the quality of leather, 
leather industry, as it leads to separation of the 
papillary section from the reticular section of the 
dermis, which affects the quality of the final products 
in leather industry. 

Cuticle i.e. flakes on the fibers from the root of the 
tail are significantly smaller, of finer structure, (Fig. 6), 
especially in regard to the ones from the rump. 

Infundibular appearance of circularly arranged 
cornified flakes may be seen on the edges of the wool 
fibers (Fig. 7). 





Fig. 5 The skin of the root of the tail (Dubska pramenka); 
the white arrow indicates dermis; the blue arrow indicates 
hair follicle; the green arrow indicates sebaceous gland; 
HXE; x 200. 
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Fig. 6 Native preparation; x 400; the root of the tail. 





Fig. 7 Native preparation; immersion; distribution of 
flakes; the root of the tail. 


This appearance is a characteristic of a finer 
high-quality wool, which contributes to high-quality 
classification bearing in mind the quality of the final 


products of the textile industry. 
4. Conclusions 


Description of histological preparations of the 
showed the 


distinctive features: thin epidermis with developed 


Dubska pramenka skin following 


dermis and hypodermis—the shoulder samples. 
Compacted collagen and elastic fibers are present in 
the papillary layer of dermis of the skin from the rump, 
stratum corneum and stratum disjunctum of epidermis 
of the skin from the root of the tail are clearly visible 
unlike in other layers. The studied samples of the 
wool fibers show certain differences in the form and 


position of the cornified flakes—cuticle depending on 


the region of the body. The cuticle of the root of the 
tail was of a finer structure; cornified flakes on the 
rump showed the features of rough wool fibers, and 
the “transitional form of the cuticle” was dominant on 
the shoulder. Distinctive features of the skin structure 
and appearance of the cuticle of the wool fibers of 
Dubska pramenka may be used as comparative 
advantages in the utilization of new raw materials in 
different industry branches and may directly influence 
improvement of production in livestock farming. 
Considering that skin has multiple functions in an 
organism and together with cuticle, it accentuates its 
protective role, it directs the further use of wool as an 


important ecological raw material. 
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Complex Use of the Water Spring Istog 
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Abstract: While living in the century of crisis (of energy and water), more focus should be given on renewable energy. Since 
Kosovo is more limited in hydro-energy resources than the neighboring countries, it is essential to study them and put to use more 
efficiently. Subject of this study is Istog spring as it is the most important spring of Mokra Gora together with Vrella and White 
Drino. The spring is analyzed as a complex resource on water economy providing: fish (trout), potable water, water for irrigation and 
hydro-energy (currently not in function). The focus of the study remains the hydro-energetic component, not just a revitalization 
process but as an upgrade of the existing facility to increase the capacity by 2-3 times. This study is based on information selection 
and processing regarding detailed technical and economic analysis providing a method for other springs that will be studied in the 
future. Two technical solutions to the problem are provided as the best economical and technical solution. Solution 1 has one level 
and a calculated flow of 7.5 m?/s and consists generally of renovation works on existing facilities; Solution 2 plans the construction 
of another level 3 m lower than the existing one, increasing so both the hydro-power capacity and implementation cost for the project. 
Both solutions provide little to no impact on the spring main attributes. The valley ecology will not be influenced because the water 
after HPP (Hydro Power Plant) will be flowing in its own bed. After this study, Istog spring will be more attractive to donors and will 
contribute in improvement of energetic structure in Kosovo that for the time being is poor in hydro component. 


Key words: Water, hydro-energy, environment engineering, complex use and sustainable management, HPP (Hydro Power Plant). 


self-effective | management services, providing 


1. Introduction—Current State of the Spring 
and Compositions on the “Istog” Spring 
Hydrology System the 


sources systems involves making decisions. Water 


unbiased and reliable public services supply that meet 


standards. Developing and managing water 


1.1 Current State of the Spring and Possible Changes 


for Rational U. sources systems typically provide a variety of 
or Rational Use 


economic, environmental and ecological services [1]. 


1.1.1 Principals, Obstacles of the Hydropower 
System 

Energy, water sectors and other public services 
suffer from long-term neglect in management, 
maintenance and investment in the physical capital. 
of the 


approaching the end of their shelf-life. They are all 


The most important parts system are 
under the theoretical capacity and there are large gaps 
in water supply and particularly in the hydropower 
system. 

based on the 


establishment principle as a stable suppler of 


The hydropower system is 
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Sciences, main research fields: environmental engineering and 
managing water resources. 


Consumers have not accepted yet the responsibility 
for payment of public services they use and that poses 
an obstacle to making progress towards the free 
market. 

De The Istog spring 1s a facility comprising of a 
few elements that have undergone several changes in 
the course of its use (Figs. 1 and 2). There have been 
changes made throughout its structure until reaching 
its current state. The schemes of the flow intake are 
performed at two levels (Fig. 2). 

The inflow composition is an arched weir of 
diameter 40.00 m, length 53.00 m, height 8.00 m. 

The catchment facility in the Istog spring is in good 
condition, both the and the 


cumulative water amount (Fig. 3). It is difficult to 


integral structure 
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515 m 
am 






Fig. 1 Out of function in the 515 m quota. 


ss 


Intake facility in use 
507.2 mM 


Fig. 2 The scheme in use in 501.2 m quota. 


Fig. 3 Photo of dam in catchment facility. 


make any changes in this facility. These changes will 
bring about the deterioration of the work regime of the 
spring itself. 

The solution required for the catchment of the water 
amount overflow is required at the lowest level, 
establishing another water catchment facility, to 
collect overflow water from the main facility as shown 


in Fig. 3. 


Old intake facility 
out of function 


1.1.2 Possible Changes on Inflow Composition 

There is a small difference between quota of the 
inflow facility and the height of the water outlet from 
the spring, which stipulates the elevation of the weir 
height, leading to the flooding of the spring and the 
deterioration of its own regime. 

In order to use the non-accessible water from this 
composition, an inflow facility should be build at a 
level slightly below the current facility (see Fig. 4 for 
two-level exploitation). 

1.2.1 Kosovo Energy System 

There is only one supplier of electricity services in 
Kosovo, KEK (Kosovo’s Energy Corporation). The 
main priorities are emergency repairs and maintenance 
to support as much as possible the power generation 
system yearly. Thus, the basis of the energy system in 
Kosovo 1s the thermo-element. 

In terms of this, it is the best opportunity to raise 
the issue of non-usage of water power in Kosovo even 
though, in comparison to e.g. Albania, water resources 


and consequently hydropower reserves are smaller. 
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Intake facility in use 


Ta. 501.2m 
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New intake facility 
438m 


Fig. 4 The new scheme that would use all the amount of water.1.2 Current State of HPP (Hydro Power Plant) and Its 


Constituent Facilities. 


Existing intake facility in the spring Channel 1 






Existing intake facility 


Fig.5 Current hydropower scheme. 


1.2.2 Current State of HPP and Its Constituent 
Facilities 

The current scheme has been built in two levels 
(Fig. 5). The HPP building is located at 474 m quota 
where the use is made at the first level with 10.00 m 
incidence. 

The current use scheme comprises channel 1, 
flowing water from catchment facility to the P.B. 
(Pressure Basin) 1 using the incidence of 380 m 
longitude. 

Catchment facility in channel 2 to the P.B. 2 using 
the incidence of 10.00 m, 135 m longitude [2]. The 
water from the two pressure basins flows to HPP 
through a pipeline. 

1.2.3 Channel of Water Inflow to the P.B. 1 

This channel performs a double function. It is 


used in agriculture between June and September and 


Existing channel 2 


Pressure Basin 1 


asin 2 





Pressure 


for hydropower the rest of the months throughout the 
year. The water amount in the channel is 2.00 m’/sec. 
The channel was totally reconstructed, bringing 1t 
to the parameters required by the two sectors (Figs. 6 
and 7). 
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Fig. 6 Photo of water transport channel. 
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Fig. 8 Photo of hydropower building. 

1.2.4 HPP Building 

From the physical standpoint, the HPP building is a 
solid structure (see Fig. 8). Considering also its 
current state of high parameters and with an internal 
architecture, capable of withstanding the further use 
and starting from the HPP building—the optimal 
facility parameters on the energy and economic basis 
will be retrieved. 

1.2.5 Pipeline Network of Water Intake to the HPP 

The entire pipeline systems are completely 
non-operational and depreciated, and they could not 
be made operational (Fig. 9). These pipelines will be 
replaced with the investments to be made. 

1.2.6 Old Turbine Type 

There are two turbines in the HPP building (Figs. 
10 and 11). They are non-operational due to serious 


depreciation and lack of maintenance. 
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Fig. 11 Properties of the existing turbine. 


These turbines are physically out of work and it 
would be impossible making them again operational 
or that could be achieved at very high costs and very 
low unjustified efficiency in the conditions of 
nowadays energy market. 

1.2.7 Channel of Water Outflow from HPP 

The channel is located on the left side of the HPP 
building and it is not concrete lined, although it is able 


to withstand the amount of water used. 
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1.2.8. Channel of Water Inflow to the Fish Farms 

Upon water outflow from HPP, the water, with an 
intake composition, flowed in through the concrete 
lined channel to the fish farm reserve (Fig. 12). 

1.2.9 Istog River Valley 

Following the HPP building and the intake water 
composition for the fish reserve, the water 
permeates quietly through the valley of the Istog river 
(Fig. 13). 

Generally, there are no erosion problems on river 
banks, which are mostly protected by trees and 
protective measures against the erosion occurrence. 
could be 
downstream where frequent flooding occurs. 

1.2.10 Possible Changes of the Facility 


In order to utilize the maximum amount of entire 


Such phenomena observed in the 


water flowing into the spring, after meeting the basic 
needs of water customers or users from this spring, 
various solutions on energy utilization of surplus 
water have been examined. 

Carrying out the necessary comparative analysis, 
the most effective solutions for the effective 
maximum use were brought to attention. 

The proposed new one-step scheme is shown in Fig. 
14. 

Two-scale scheme 2 (Fig. 15) by changing the 
current scheme, the first scale at 501.2 m quota with 


the current channel and the second scale at 498 m 


Channel 1 


Existing intake facility in the spring 





Existing intake facility 
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Fig. 14 Scheme with one degree of exploitation. 


Existing channel 2 


quota with steel pipeline with 315 m length is shown 
in Fig. 15. 





Fig. 12 Reserve of trout in Istog. 





Fig. 13 Istog river valley. 
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Fig. 15 Scheme with two degree of exploitation. 
2. Hydrological Study on the “Istog” Spring 
2.1 Hydrology of the Istog Spring 


The hydrological study aims at recognizing and 
evaluating the overall water potential of the country 
and the laws of this potential distribution in space and 
time. The utilization of the hydrographic water 
network for the needs of economic and cultural 
development must be preceded in order to recognize 
this value, without allowing the abundant waters of 
streams and rivers, being a gift of nature, to flow in 
vain, without utilizing their usefulness for hydropower, 
irrigation, water supply of population and industry, 
recreation etc.. 

The aim of this paper, relying on the processing and 
homogenization of a multi-annual comprehensive 
hydrometric material for the period 1938-1992, was to 
assess both the watercourse of the Istog underground 
spring and its relation to the physical-geographic 
factors of the catchment basin. The watercourse is the 
most important element of a water resource. It also 
determines, better than any other element, the ability 
of the resource to be used for different purposes. 

Waters of the hydrographic network of the Istog 
spring are widely used for drinking water supply to 
population, irrigation of agricultural lands, industry, 
energy, fish farming, etc.. 

2.1.1 
Catchment Basin of the “Istog” Spring 


Physical-Geographic Conditions of the 
The Istog spring is one of the main underground 


springs of the Kosovo hydrographic network. The 


Existing intake facility 


—_ i 
New intake facility F = = = = 


Pressure Basin 1 







New pressure basin 3 









hydrological regime of the spring is closely related to 
physical-geographic characteristics of the catchment 
basin. 

The main physical-geographic characteristics of the 
catchment basin, affecting the hydrological regime of 
this spring are as: 

e Geographic expansion to the catchment basin; 

e Relief; 

e Hydrography; 

e Geology; 

e Pedology; 

e Tectonics; 

e Vegetation cover; 

e Geographic distribution of climate elements. 

2.1.2 Geographic Expansion to the Catchment 
Basin 

The catchment basm that forms the Istog Basin 1s 
located on the northern part of Dukagjini plateau, 
between the parallels 42°40’ and 42°70’ of the 
northern latitude and the 30°21’ and 30°20’ meridian, 
42° eastern latitude [3]. It includes approximately 1/9 
of the entire surface of the plane. The northern part of 
Istog river basin lies in Accursed Mountains and a 
good part at Dry Mountains. Its altitude reaches 1,758 
m at the peak grape, followed by Dry Mountains in 
the north-western part, with the highest peak of 2,154 
m above the sea level. It lies entirely in the territory of 
Kosovo, with the general flow course from north-west 
to south. Located in the south of the “Mokra Gora” 
Mountain (the Wet Mountain), the Iber river lies in the 


north over the mountains, in the underground spring 
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Vrelle. The basin is 53% mountainous and 47% field 
area. The surface of the catchment basin is 440 km” 
and the river basin is 216 km?. The basin has very 
convenient connections to the second and third 
category routes, located in the southern and central 
part of the basin. 

2.1.3 Relief 

The configuration features of the catchment basin 
of the Istog spring are determined by the geographical 
expansion of the basin area. The altitude of the spring 
in 520 m quota above the sea level shows that this 
spring collects waters of a territory with an 
emphasized mountainous character. The highest point 
is about 1,500 m above the sea level. The expansion 
of the basin creating the Istog underground spring is 
76 km’, and altimetric position over 1,500 m. Its area 
is 52 km? at 1,000-1,500 m quota about 22 km’, 
under 1,000 m with 2 km? area, karstic part of the 
catchment area of 6 km’, feeding the Istog spring [4]. 

2.1.4 Hydrography 

The hydro-graphic area of the Istog spring 
catchment basin consists of a completely underground 
network. This area is characterized by small springs 
flowing to its main effluent, which is the Istog spring 
itself and consisting of two adjacent springs. In 
various periods, there are fluctuations of water 
quantities in these springs. During the period of inflow 
measurements, Q = (1.19-4.9) m'/sec are observed, 
the minimum inflow rate is recorded between 
September to October 1933 (1.6-1.8) m'/sec which 
continued for about 50 days and the maximum inflow 
in 1936 is about 6.6 m'/sec [5]. The annual volume 
flowing into Istog is (75-126) x 106 m°. The use of the 
spring is complex. It is used for drinking water supply, 
agriculture and trout reserve. It was also used for 
hydro-electricity in the past. 

In regards to the water used for HPP, there is a 
possibility in taking more water and at higher levels. 

The Istog spring is the first spring in the Dukagjini 
hydrographic network, joining other effluents which 


are generally underground springs of the Vrella spring, 


the White Drino spring, Bistrica of Peja, Bistrica of 
Decani, creating the White Drino river, which 
percolates to Albania. 

This network is considered as the richest one with 
groundwater basins. The comparison of hydrological 
conditions of the hydrographic network shows that 
they have made the accumulation and rational use of 
land protection from 


water resources possible, 


flooding, water systematization in terms of 
agricultural intensification, etc.. 

2.1.5 Assessment of Hydrology in the Catchment 
Basin Forming “Istog” Underground Basin 

As a result of climatic factors, the regime of flow 
distribution throughout the year in the Istog spring is 
largely formed by rainfalls, soil surface factors, 
mainly by the relief of geological formation. 

The first observations carried out by the department 
of Peja water measurements belong to the years 
1933-1940. Measurements were made by an engineer 
called Burovic, whose values are reliable because the 
measurements were made with overflow facility. 

Figs. 1 and 2 show annual flow course of the years 
1933 to 1940, 


measurements. 


and 27 years for recorded 

Fig. 3 shows the evolution of the average flow of 27 
years. 

Maximum and minimum values were recorded 
from the observations made for the years 1933-1940 
(Table 1) and 1957-1959 (Fig. 16). 

Maximum annual inflow Qmax = 4.0 m'/sec [6]. 

2.1.6 Geology 

From the previous studies, it was concluded that 
this pond was built entirely of thick plaque limestone. 
Karst springs were formed as a result of geological 
formation. This composition has conditioned a radial 
extension of the small karst pits and troughs of 1,500 
m thickness, on which a flavio-dolomite broad plane 
is formed. 

Two corrugated directions were combined in this 
section: the direction of characteristic abbreviations 


table on Accursed Mountains and karstic rocks direction, 
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Table 1 Average monthly inflow 1933-1940. 
Monthly avg. inflow 1933-1940 avg. 3.09 m°/sec 
Month I II HI IV V VI VII VU IX X XI XII 
Inflow (m°/sec) 2.72 275 2.96 3.19 É a 3.53 do AZ 2.84 2.88 3.05 3.09 
Ag 4 
® ee ee S 
V 3 
(e) 
E 2 
= 
5 1 
Sam 
a 0 
| II || IV V VI VII VIII IX X xI XII 
—e—m3/sek 2.72 2.75 | 2.96 3.19 3.57 3.53 3.3 3.12 2.84 2.88 3.05 3.09 
Fig. 16 Average monthly flows 1933-1940. 
which is known as dynamic management, the erosion process, as due to the superficial 


characteristic for the Paleozoic formations. In general, 
they dominate the limestone rocks composed of 
MgCO; magnesium, Ca calcium and dolomite. 

This block represents the highest and the most 
magnificient part of “Dukagjini”, the Kosovo Alps, 
situated in the center of them. 

2.1.7 Pedology 

The surface of the catchment basin of the Istog 
spring, for its own meridian extent and rugged relief, 
in terms of pedological aspect, is formed of different 
types of soil. 

Pasture and carbonated-humus lands are located on 
the altitude of 1,000-1,200 m above sea level. These 
lands, mechanically placed under the pasture alpine 
and subalpine vegetation, represent the upland surface 
source. These soils occupy a small area of the total 
territory of the catchment basin of this spring. Brown 
forest lands, stretching at 700-1,000 m height above 
sea level, that are mechanically placed under beech, 
chestnut, Mediterranean fir vegetation etc. represent 
about 18% of the forested area. 

These of high 
hydro-capacity. These lands are distinguished by low 


soils are permeable and 


content of humus. For the most part, the soil cover or 


very 
Consequently, these lands are more easily subject to 


the humus layer is shallow or absent. 


mechanical action, they detach more easily from the 
surface of the catchment basin, therefore, large 
infiltration occurs due to this phenomenon. As a result 
of the pedologic construction, the catchment basin that 
forms the Istog spring is generally easily subject to the 
erosion process. 

2.1.8 Tectonics 

The Dukagjini Plain (Rrafshii Dukagjinit) has 
gained the shape of its territory in the newer era of 
geological formation. Its forests started to be formed 
in the Diogenes corrugation towards east-west and 
The 


predisposition has defined the basic contour elements. 


north-south-west alignment. tectonic 
Radial movement in this region is reflected in 
orogenies of the mountainous territory to the west 
and south. Small and large slides are observed in this 
date the 


tectonic-geomorphology shape has seen many major 


area. Since its creation to 
modifications and they have lost the shape of the 
previous initial relief. 

Parts of these mountains are not as rugged and 
harsh as before. Now they are more rounded and 
increasingly losing their previous shapes. 

2.1.9 Vegetation Cover 

Another characteristic is also the scarce vegetation 


Cover. 
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By vegetation cover, the surface of the catchment 
basin forming the Istog spring is divided into different 
phytoclimatic areas. 

Phytoclimatic area of subalpine pasture is situated 
in the catchment basin territory and located at 1,000 m 
altitude above the sea level. 

Shrub area lies at lower altitudes, or in the lower 
area of the territory. 

Beech, oak, phytoclimatic areas lie at 900 m 
altitude, and the chestnut, white pine area lie at 
700-900 m above the sea level. 

This area covered by vegetation continues merging 
with the surface of the catchment basin of Vrella. 
These two underground basins are the main effluents 
forming the underground basin Lubozhde. Therefore, 
most of the trees are deciduous and protect land 
slightly from the action of meteorological agents. 
Forest massifs are mostly poor, they are rare, not 
uniformly distributed across the surface, and 
consequently, they are not able to collect large 
quantities of water and make the evaporation process 
intensive. 

2.1.10 Geographic Expansion of Climatic Elements 

2.1.10.1 Climatic Position 

Due to its mountain extent, the surface of the 
catchment basin falls under the humid sub-tropical 
belt. Significant climate changes have been observed, 
expressed in the gradual transition to the average 
continental climate zone of the Central Europe. The 
main elements conditioning the climate catchment 
regime as precipitation, air temperature, humidity, etc. 
are characterized by an annual course more or less 
designed, depending on the height of land above the 
sea level, the geographical scope of the basin from the 
intensity of atmospheric circulation. The area is 
characterized by cold winters and long duration. The 
snow there is an ordinary occurrence, while the height 
of snow layer reaches up to 2 m. 

2.1.10.2 Air and Water Temperature 

Air temperature is one of the mam factors 


characterizing the thermal regime of the basin. 


Analyzing the dependence of the thermal regime with 
physical-geographic characteristics of the 
hydrographic basin territory, a conclusion that in 
general, there is a connection to some extent between 
the altitude of the country and the height of air 
temperature is reached, however, the character is 
regional in nature. Average lowest annual 
temperatures observed are 11.3 °C. The main factor of 
temperature change in this area is altitude and 
geographical position. It turns out that the coldest 
month is January and the warmest is June-July. 

Average annual water temperature is 7-9 °C. 

2.1.10.3 Rainfalls 

They are main hydro-meteorological characteristic 
and their distribution in space and time greatly 
influence the flow regime, because it directly affects 
filling up water reserves on the land surface and 
underground. The main features of the pluviometric 
regime of a geographic region can be recognized by 
the average annual values. The part where the Istog 
catchment basin lies is included in the area 
characterized by the largest average annual amount of 
rainfalls, where annual rainfall reach up to 1,500 
mm/year. The largest amount of precipitation falls in 
December and the lowest in June. 

Precipitations in the form of snow fall mainly 
during winter months, and their sustainability depends 
on local weather conditions and air temperature. 

2.1.10.4 Solar Radiation 

Data analysis shows that the intensity of the total 
radiation for unclouded weather conditions has regular 
annual performance. The maximum intensity value is 
in July, and the minimum value in January. Total 
maximum radiation along the territory is 24.79 


kWh/m’. Minimum falls down to 24.79 kWh/m”. 


3. Environmental Impact 


One of the contemporary issues in the construction 


of engineering works and in particular of 
hydro-technical works is the environmental impact 


study. 
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different 


methodologies developed to assess the environmental 


In western practice, there are 
work impact, amongst which is the NNC (National 
Research Center—Italy) under which the environment 
is studied by components and factors according to 


Section 3.1. 
3.1 Natural Constituents 


3.1.1 Physical Constituents 

The physical constituents are surface waters, 
ground waters, soils, climate, lithology, geology, air 
and noises. 

3.1.2 Biological Constituents 

The biological constituents are mammals, fish, 
amphibians, reptiles, birds, insects, microorganisms, 


herbal vegetation, shrubs, natural vegetation and crops. 
3.2 Socio-economic Constituents 


The socio-economic constituents are population, 
exodus, culture, archeology, residential zone, health, 
safety, industry, trade, services, tourism and 
agriculture. 

This 1s a current and present methodology in 
Europe and widely, however it is less known and not 
used at all in Kosovo and Albania. 

Certainly, weirs and reservoirs hydro-technical 
actions make a contribution to the environment, 
however, derivation hydropower actions also have an 
impact, which should be assessed both 1n positive and 
negative terms. 

The essence lies in the | assessment 
causality-consequence of the interaction 
work-respective territory. 

As seen from Section 3, the main components are 
component of the natural environment (includes 
physical and biological component) and 
socio-economic component. 

These components have their own factors, which 
should be identified and assessed before and after the 
implementation of the project. 


What is important in this methodology is the use of 


assessment “weight attribute” criteria factor, namely, 
the impact of the project on each factor is assessed not 
only for quality, but also for quantity. Thus, each 
factor is given appropriate relevance in relation to the 
others, reaching finally an “IQUAM (Final Weighted 
Environmental Quality Index)”. 

The main procedural stages of an environmental 
compatibility study shall be: 

(1) Preliminary researches; 

(2) Description and identification of actions and 
factors; 

(3) Weighing attribution; 

(4) Impact assessment; 

(5) Data formatting; 

(6) Final assessment. 

Meanwhile, in order to reach to IQUAM, the genric 
environmental structure is developed at four levels as 
shown in Table 2. 

In order to assist the methodical application, the 
ICOLD (Italian National Commission on Large Dams) 
has also made the list of designing actions and use of 
compositions. This list shall be used as a standard in 


the EIS (Environmental Impact Study). 
3.3 Water Use 


The water use are irrigation, energy, drinking 
water, regulation of rivers, industrial use, sailing, use 
against fire, fishing, sport use and environmental 


various use. 
3.4 Physical Factors 


The physical factors are weir presence, reservoir, 
water diversion, construction site, dikes, buildings, 
deforestation, closure, channels and pipeline, water 
outflows, inflow and diversion actions, transmission 


lines and unconventional materials. 


Table 2 Generic environmental structure levels. 


At level zero General environment 


At level I X, Y, Z systems. 
At level II Components for every system 
At level HI Factors for each component 
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3.5 Areas of Interest 


The areas of interest are flooded areas, lake banks, 
water level oscillation in the lake, upper part of river, 
lower part of river, irrigation channels, ground waters, 
sea coast and lower part of river connected to the 


basin. 
3.6 Corrective Actions 


The corrective actions are phology repopulation , 
guaranteed river flow, tourism development, water 
level control, infrastructure, reforestation, erosion 
control, river deepening, weir, outflow, auxiliary weir, 
compensation reservoir, barriers against floating 


bodies, 


improvement industries, transfer of population and 


basin protection, water treatment, 


area of protection from avio fauna. 
3.7 Normative Predictions 


The normative predictions are national laws, 
regional laws, urban regulations, contributions and 


taxes and immovable property values. 
3.8 Interrelation Factors 


The interrelation factors are: inter-regional 
reservoirs, inter-municipal reservoirs, benefits for 
other regions and benefits for other municipalities. 

The quality impact assessment is also performed in 
the order shown in Table 3. 

Authors emphasizes that the application of this 
methodology (in particular the section of relevance 
assessment—weighing of each factor or component 
that is based on “pair” factors confrontation under the 
formulated technique by SAAT (Sanitation Approach 
Advise and Training) is not simple and certainly, its 
depth depends also on the stages of the project. 

In this case study, several possible estimations have 
been made, especially qualitatively, without analyzing 
the quantitative part of the environmental impact. 
Certainly, this action also requires fulfillment in 
compliance with the proposed methodology, however 


this will be implemented in further stages of the 


Table 3 Quantity impact assessment 


Type of impact Symbol 

Impact Positive E 
Negative - 

Effect Irreversible P 
Reversible K 

Duration Continuous V 
Temporary P 

Area of impact Zonal Z 


study—design on the complex use of this water 


source. 


4. Energy Development Strategy 


The strategy pursued in the public service sector 
aims to replace the periodic cycle in emergency aid to 
a gradual process of a sustainable development. The 
main goal in the public service sector is the 
establishment of sustainable providers of public 
services with effective self-management in defense 
and the public interest services by providing impartial 
and credible public procurement, meeting proper health, 
safety, environmental services standards with reasonable 
prices for all customers. An important and quantifiable 
criterion for plan is the economic benefits and costs a 
plan would entail when it implemented [7]. 

There is only one provider of electrical services 
(KEK) in Kosovo. The top priorities are urgent repairs 
and maintenance to support as much as possible the 
power generation system for the winter period. There 
have been extensive studies on long-term capital and 
operational opportunities and future needs. Detailed 
reviews have been made on the overall state of the 
system in all aspects. 

Planning has triggered discussions and mid-term 
and long-term projects on the initiation of a project—a 


comprehensive study covering a long-term period. 
4.1 HPP Location 


The composition actions of the hydropower node 
are located on the left side of the Istog river valley 
(Fig. 17). The inflow composition is 501.2 m and at 
498 m quote, the channel is directed to 501.2-500.2 
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Fig. 17 Location of the hydropower plant. 


m and the pipeline of water delivery to P.B. 3 basin is 
at 498-497 m quota. The HPP building is situated at 
an altitude of 474 m and is about 4.00 m above the 


river valley. 


5. Environmental Impact, Assessments and 
Counteraction Measures 


5.1 Structures and Roads 


There were no river deviations during the works 
which do not affect the valley ecology. 

Roads are not affected in this project. During 
construction and reconstruction, the facilities of other 
sectors, agriculture, water supply and the fish (trout) 
farming plant are not affected at all, which inflow 
water from the HPP water outflow channel. Making 
the HPP again operational would have a positive 
effect on the Istog area and beyond because the source 
is simultaneously fully utilized, electricity is produced 


and a significant number of workers are employed. 
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5.2 Visual Effect 


The visual effect is not affected (Fig. 18) and the 
reason is that there are many constituent compositions 
that will be reconstructed with the exception of the 
second instance that will be a new composition. The 
solution of pipeline water inflow to the pressure basin 
has been recognized in order to maintain the visual 
effect. 


5.3 Flora and Fauna 


Situated in a mountainous area, flora and fauna are 
the most developed elements both throughout the river 


valley and forest slopes. 
5.4 Land Use 


Lands in the Istog river valley are in good physical 
conditions, cultivated in all surfaces, systemized, 
equipped with all necessary irrigation and drainage 


systems. 


526 Complex Use of the Water Spring Istog 





Fig. 18 Visual effect. 
5.5 Population and Employment 


The town of Istog with its surroundings has a 
relatively small population of 10,000 inhabitants. 
Employment and labor market is limited because the 
town of Istog has a limited physical position. 

Restoring the HPP, fishing and the development of 
agriculture, these sectors reduce the unemployment 


level and raise the living standards. 


6. Conclusions and Recommendations 


From the conducted study, it can be concluded as: 

e The results achieved from this study support the 
general trend in the world to reach a better 
understanding towards better and complex use of 
renewable resources. Water resources management is 
an iterative process of integrated decision-making 
regarding the uses and modifications of waters and 
related lands within a geographic region [8]; 

e Although Kosovo, compared to the neighbor 
countries, is poor in water resources (hence 
hydropower resources), it is an urgent task starting 
and conducting studies as soon as possible on all 
streams, progressively drafting a master plan in now 
new conditions, treating water as a vital element in 
order to avoid emerging crisis 1n the near future; 


e With the presented solutions (without affecting 


their components of water economy: on water supply 
of the population, agriculture, fish economy, valley 
ecology), the Istog source may provide an average 
annual electricity production of nearly 7 million kWh 
with a gross income of about 300,000 euro per year; 

e The implementation of this study is extremely 
easily realized for a very short term; 

e It is recommended to colleagues who are in 
various stages of their doctorate to be involved in 
similar studies for other water sources, specifically 
contributing to gradually complete an overview of 
water and hydropower opportunities of the state of 
Kosovo. To this regard, the department, faculty in 
cooperation with relevant local government ministries 
have a real opportunity to contribute to this nature; 

e Authors of this study is aware that his work has 
flaws but it is a step 1n the right direction in paving the 
road for conducting concrete studies and serving 
Kosovo’s economy, creating a foundation where 
following works will be more profoundly performed 


from a scientific standpoint. 
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Estimation of Soil Hydrodynamic Parameters Related to 
Agricultural Practices —Case of the Tougou 


Experimental Site (Burkina Faso) 
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Abstract: This study aims to estimate the hydrodynamic properties of soils under various agricultural practices in the Tougou 
catchment in northern Burkina Faso. The methodology adopted is based on the determination of the unsaturated hydraulic conductivity 
and capillary sorptivity close to saturation. This method relies on the measurement of the transient infiltration flux at the soil surface 
with imposed hydraulic head varying from -60 to -20 mm. These tests are carried out on control, stony line, half-moon and zai plots. 
The results show a difference in hydrodynamic parameters according to the agricultural practices. The unsaturated hydraulic 
conductivity is 33.1 cm/h, 13.1 cm/h, 20.3 cm/h and 4.0 cm/h for zai, control, stony line and half-moon plots respectively. The 
unsaturated hydraulic conductivity is 33.1 cm/h, 13.1 cm/h, 20.3 cm/h and 4.0 cm/h for zai, control, stony line and half-moon plots 
respectively. The pores participating to water transfer also differ. The mean size of drainable pores is 43.7, 56.2, 22.3 and 87.2 um on 
control, stony line, half-moon and zai plots respectively. 


Key words: Agricultural practice, unsaturated hydraulic conductivity, sorptivity, soil, Tougou, Burkina Faso. 


1. Introduction of rainfall (splash effect) and runoff. It thus appears at 


landscape scale an imbalance in the spatial distribution 
The Sahel of Burkina Faso, like all the Sahelian = 
of water, accentuating the constraining nature of this 
regions, undergoes for several decades various 
resource. To this, added anthropic pressures and the 
phenomena of natural and anthropogenic degradation 
development of livestock farming thus increasing the 
of ecosystems. This degradation, accentuated by a o 
l l vulnerability of the environment to drought [1]. 
decrease in rainfall over the last decades, causes a , , 
, , , The degradation of natural resources experienced by 
modification of the processes and mechanisms that , a o 
Sahelian countries in recent years limit their capacity 
govern these ecological units, which ultimately leads to 
for endogenous development. To increase their 
a modification of the state and aspect of these o , , 
, , , resilience to this, Sahelian farmers in Burkina Faso 
ecosystems. The various anthropogenic actions , , 
, i , , have adopted several agricultural practices for the 
combined with the rainfall deficit, cause a reduction of ee 
or water mobilization at the farm scale [2-7]. These 
the vegetal cover, or even locally its disappearance. 
practices, which appeared between 1980 and 1985 and 
Indeed, in areas where vegetation has disappeared, it oo, 
are very demanding in terms of working time and 
has developed on soil surface an indurated, continuous 
physical effort, nevertheless constitute innovative 
and impermeable layer (superficial crust). This layer ee 
aa alternatives for the development of a resilient 
set up is controlled by the action of the kinetic energy l , , , 
agriculture in the Sahel. They consist of applying 
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fields: hydrology and soil physics. line and using of short cycle varieties. These 
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techniques have allowed not only the restoration of 
many uncultivated lands [8], but also the improvement 
of cereal yields [9] and the increase of the groundwater 
level [10]. Fatondji, D., et al [8] reported that with the 
zai technique 100,000 to 300,000 ha of degraded land 
were rehabilitated in Burkina Faso while Roose, E., et 
al. [9] showed in the same country that under zai, 
sorghum yields are multiplied by about 1.6. These yield 
increases result from the combined effect of water 
retention, composting and supplemental manure 
fertilization. 

In such context, knowledge of soil water transfer 
processes 1s a key priority. Thus, the present study is 
implemented on an experimental site in Tougou 
located in the north of Burkina Faso. It aims at 
estimating soil hydrodynamic properties related to 


several types of agricultural practices. 
2. Material and Methods 
2.1. Study Area 


The study area is the Tougou catchment (between 


Latitudes of 13°40’ N, and longitudes of 2°13’ W), a 
surface sub-catchment of Nakambe river in the 
north-east of Burkina Faso. It has a surface area of 37 
km? (Fig. 1). The climate is semi-arid with a mean 
annual rainfall varying from 400 to 650 mm. 
Temperatures range between 18 and 40 °C. The dry 
season extends from October to May and the rainy 
season from June to September with peak precipitation 
generally recorded in July or August. Soils are 
cultivated or denuded (degraded) and the vegetation 
consists of savanna, shrub and grassy steppe. 

The experimental design (Fig. 2) consists of two 
separate Fisher blocks A and B, the distribution of 
which is random. Each block consists of four plots of 
200 m? (20 m long by 10 m wide) each receiving a 
specific treatment: direct seeding as a control plot (TO), 
a stony line plot (T1), a half-moon plot (T2) and a zai 
plot (T3). Zai (plot T3) is a complex system of 
restoring the productivity of degraded soils, a particular 
type of an open hole (20-40 cm in diameter and 10-15 


cm deep), concentrating runoff and organic matter. The 
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Fig. 1 Location of Tougou surface catchment. 
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Fig. 2 Design of agricultural practices plots in Tougou experimental site. 


spacing between the zai holes is 60 cm x 60 cm. The 
stony lines (plot T1) are constructed as a hydraulic 
barrier of polyform stones anchored 15 cm deep and 
located 15 m from plot upstream limit, parallel to the 
contour lines to dissipate surface runoff energy. 
Half-moons (plot T2) are semicircular holes of 4 m in 
diameter and 15 to 25 cm deep, and are excavated 
perpendicular to water flow direction. The spacing 
between half-moons is 2 m. The control plot (TO) is a 
direct seeding, in rows, as practiced in the locality. The 
spacing between seedlings on the control, stony-line 
and half-moons plots is 50 cm x 50 cm. An organic 
amendment (cow dung) of 5t/ha is applied to each plot 
prior to seeding and a microdose application of NPK 
(Nitrogen, Phosphorus, Potassium) fertilizer (2-4g per 
hole) after crop emergence and urea (1 g per hole) at 
the run [11, 12]. Plots are grown with a 70 days cycle 
millet variety (Kiipalla). The choice of this variety is 
justified by its precocity, its resistance to drought and 


its wide adoption by local farmers. 
2.2 Field Data Collection 


Measurements of the unsaturated hydraulic 
conductivity are implemented using disc infiltrometer 
on all eight plots of the experimental site. On each plot, 
the infiltrometry tests are carried out at suctions 20, 40 
and 60 mm with a repetition of 3 measurements per 


suction and the average value was recorded. Before 


starting and ending each test, soil samples are analyzed 


in order to estimate the initial and final water contents. 
2.3 Data Processing 


2.3.1 Unsaturated Hydraulic Conductivity and 
Capillary Sorptivity 
The three-dimensional cumulative infiltration per 
unit of area, I (mm) for the entire time range, can be 
expressed by Eq. (1) [13]. 
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where Rp (m) is the radius of the disc; Oo and 0, are the 
final and initial volumetric water content (m°-m°>), 
respectively; So is the sorptivity (m-s””) for 69; and y is 
the proportionality constant, the value of which can be 
approximated to 0.75 [14]; Ko and K, are the soil 
hydraulic conductivity values (m-s!) corresponding to 
Oo and 0n, respectively; and B is a shape constant that 
commonly takes an average value of 0.6 [14]. In spite 


of its relative complexity, Eq. (1) is valid for the entire 
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time range, from t=0 to t= œœ. However, taking 
into account that infiltrometer experiments do not 
require very long time of application, 
Haverkamp, R., et al. [13] established that, for short to 


medium time and assuming K, = 0, Eq. (1) can be 


ranges 


simplified to 
Ign=CyVt+Cot (2) 
Where 
Cy = So (3) 
2 = vS6 
C, = —— K + (4) 
: 3 Rp(6n-00) 


Using this expression, Vandervaere, J., et al. [15] 
proposed the differentiated linearization method to 
infer soil hydraulic properties using linear regressions. 
The technique consists in differentiating Eq. (2) with 
respect to the square root of time. 

dl 


dyt 


And next plotting the 


= Ci + 2C,vVt (5) 
ar 
dvt 
The Cl is the intercept and C2 is the slope of the 


corresponding regression lines. According to the 


term as a function of Vt. 


authors, the differentiated linearization technique 
allowed visual monitoring of the contact layer, when 
used, eliminating its influence on the estimates of the 
soil hydraulic properties. 

2.3.2 Mean Size of Drainable Pores 

The knowledge of the capillary sorptivity SO and of 
the hydraulic conductivity K has been used by some 
authors [16, 17] to define the capillary length Ac which 
is a macroscopic scale expressing the relative 
importance of the capillary and gravitational forces 
acting on water penetration in the soil. Its mathematical 


formulation is as Eq. (6): 
Ae = pe (6) 


Where b is 0.55, a value generally accepted for most 
soils [18]. Using the elementary laws of capillarity, 


Philip, J. R. [16] introduced the mean size drainable 


pores, Am, given by Eq. (7): 
TE 
m Pw Ih 
Where o is the surface tension coefficient of water 
(0.072 N/m at 25 °C), py the density of water (1,000 
Kg/m’) and g the acceleration of gravity (9.81 m/s). 


(7) 


By introducing in Eq. (6), the numerical values of o, py 
and g are obtained by expressing Am in um Eq. (8): 
(On a 0o)K 


ba =I z? 


(7) 


3. Results and Discussion 


3.1 Results 


3.1.1. Unsaturated Hydraulic Conductivity and 
Capillary Sorptivity 

The unsaturated hydraulic conductivity and the 
capillary sorptivity according the suction follow a 
logarithmic law and change in a decreasing manner 
when the suction increases on all the plots. As an 
illustration, the evolution on the plot T3A of these two 
parameters as a function of the suction is represented in 
Fig. 3. 

On each of the two blocks of the experimental site, 
the highest unsaturated hydraulic conductivity values 
(Table 1) are found in the zai plots and the low values 
in the control plots. With regard to the sorptivity values, 
the half-moon plots show the highest values whereas 
the zai plots show the lowest values. 

3.1.2. Mean Size of Drainable Pores 

The distribution of mean size of pore as a function of 
agricultural practices (Fig. 4) indicates that the 
maximum of pores participating to the flow is recorded 
on zai plots. This is totally in agreement with the results 
obtained for the unsaturated hydraulic conductivity. 
This is certainly due to the amendment in organic 
matter made on this practice which allows increasing 
the microbiological activity of the soil. The latter tends 
to favor the creation of a preferential flow path for 
water in the soil, hence the importance of the 


infiltration capacity of this practice. 
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Fig. 3 Distribution according to the suction of unsaturated hydraulic conductivity and capillary sorptivity on plot T3. 


Table 1 Unsaturated hydraulic conductivity (K) and capillary sorptivity capillaire (S,) recorded at a suction 40 mm. 


Plot K (mm/s) 
TOA 0.031 + 0.008 
TIA 0.048 + 0.005 
T2A 0.009 + 0.001 
T3A 0.086 + 0.012 
TOB 0.042 + 0.008 
T1B 0.065 + 0.010 
T2B 0.013 + 0.0002 
T3B 0.098 + 0.012 


So (mm/s"”) 


0.052 + 0.003° 
0.042 + 0.008 
0.067 + 0.010 
0.029 + 0.002 
0.040 + 0.005 
0.034 + 0.006 
0.044 + 0.004 
0.021 + 0.005 


a: Variation coefficient for hydraulic conductivity; b: Variation coefficient for capillairy sorptivity. 


3.2 Discussion 


The infiltrometry tests demonstrate that unsaturated 
hydraulic conductivity and sorptivity decrease with 
increasing suction. Indeed, these two parameters which 
follow a logarithmic law, vary with soil water content 
and the suction. At suctions close to saturation, the 
movement of water is controlled by the coarsest 
macroporosity. The soil has a high hydraulic 
conductivity associated with high porosity. As the soil 
becomes unsaturated, the larger pores are emptied and 
the flow occurs in smaller and smaller pores with more 


and more tortuous flow paths. These results confirm a 


pronounced differentiation of these soils according to 
the agricultural practice with the consequence of an 
accentuation of the spatial discontinuity, as well from 
surface roughness as of their hydraulic behavior. In 
addition to highlighting this diversity of surface 
conditions, the analysis of the results obtained with 
regard to the hydraulic properties of the surface reveals 
a strong variation of the infiltration capacity. The 
highest values are found on the zai plots while the 
lowest values are found on the control plots. These 
are undoubtedly attributable to the 


reorganization of land cover by the implementation of 


differences 


the agricultural practices because a modification of the 
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Fig. 4 Distribution according to agricultural practices of mean size of drainable pores. 


organization of the poral system of the superficial 
horizons leads to a reduction of the hydraulic 
conductivity close to saturation [15, 19, 20]. The 
hydraulic conductivity values estimated on the control 
plots correspond to those found on degraded soils. It 
reflects the hydrodynamic behavior of soils consisting 
generally of glacis having a superficial layer 
characterized by a low infiltration capacity due to a 
reduced porosity, which appears this laminate coating 
to a real hydraulic barrier which strongly limits the 
water inflow in the soil. This type of land cover 
generates large runoff and correlatively leads to a very 
low accumulation of water in the soil. The hydraulic 
conductivity values estimated on the zai and stony line 
are larger than those of the control plot. Both of these 
agricultural practices tend to improve the infiltration 
capacity of the soil. For zai, the destruction of the upper 
soil layer will favor the capture of runoff, its solid and 
soluble load, and tends to increase the infiltration 
capacity [21]. This leads to high hydraulic conductivity 
values resulting from the improvement of the soil 
Indeed, 


conservation and the concentration of nutrients will 


roughness. runoff management, manure 


create conditions for optimum water transfer in the soil 


through the generation of preferential flow paths (high 
termite activity). Concerning the stony line plots, 
improvements in hydrodynamic parameters are not 
very significant compared to the control plot in that the 
experimental design is in its second year of existence 
and several studies [22, 23] carried out in this region of 
Burkina Faso showed that such practice becomes 
efficient after three years. The results obtained on the 
half-moon show very low values compared to those of 
the control plot. This is certainly due to the fact that this 
agricultural practice tends to carry a lot of suspended 
matter which during the settling phase will clog the soil 
thus making the infiltration capacity very weak. The 
comparison of infiltrometry tests show that unsaturated 
hydraulic conductivity and sorptivity do not always 
change in the same way. In fact, the decrease of the 
hydraulic conductivity is much faster than that of the 
sorptivity. This behavior is due to several factors 
among which one can mention the presence of a 
vesicular porosity which reduces the capacity of 
infiltration. This vesicular porosity appears following 
the rains and is formed by trapping air during heavy 
rainfall. Indeed, some authors [24] have shown that in 


arid zone, the hydrodynamic characteristics of 
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superficial microhorizons could be affected by the 
appearance of these vesicles. This vesicular porosity is 
linked to the existence of an underlying microhorizon 
with low porosity, most often a material with low gas 
diffusivity [24]. This results in a reduced hydraulic 
conductivity of the superficial organizations. However, 
several authors have shown that in the Sahelian zone, it 
is the hydrodynamic characteristics of superficial 
microhorizons which condition, to a large extent, the 
infiltrability of the whole soil. It is therefore possible to 
consider the presence of these vesicles close to the 
surface as an index of low permeability [24] because 
the statistical analysis shows a very good relationship 
between the abundance of these pores and the runoff 
ability [25]. 

Numerous field or laboratory studies on undisturbed 
soils have shown a very significant increase in 
hydraulic conductivity close to saturation. This 
increase is generally greater than that expected from an 
exponential relationship of K (h), a relationship often 
verified m soil physics for remolded soils or with a 
particularly homogeneous structure [26]. Jarvis, N. J. 
and Messing, I. [27] have shown that, for several soils, 
this increase was particularly marked for water heads 
greater than -60 mm, i.e. when pores of equivalent 
radius greater than 0.25 mm participate to the flow. 
These authors attribute these observations to the effect 
of soil macroporosity, defined as the set of pores with 
an equivalent radius greater than 0.25 mm (thus visible 


to the naked eye) participating to water transfer. 
4. Conclusion 


This study revealed that unsaturated hydraulic 
conductivity and sorptivity decrease with increasing 
suction. The hydrodynamic behavior of the soils under 
agricultural practices reveals a strong variation of the 
infiltration capacity. The highest values are found on 
the zai plots while the lowest values are found on the 
control plots. The unsaturated hydraulic conductivity is 
33.1 cm/h, 13.1 cm/h, 20.3 cm/h and 4.0 cm/h for zai, 


control, stony line and half-moon plots respectively. 


The pores participating to water transfer also differ. 
The mean size of drainable pores is 43.7, 56.2, 22.3 and 
87.2 um on control, stony line, half-moon and zai plots 


respectively. 
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